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Abstract 

The analysis of seasonal or annual block maxima is of interest in fields such as 
hydrology, climatology or meteorology. In connection with the celebrated method 
of block maxima, we study several tests that can be used to assess whether the 
available series of maxima is identically distributed. It is assumed that block max¬ 
ima are independent but not necessarily generalized extreme value distributed. The 
asymptotic null distributions of the test statistics are investigated and the practical 
computation of approximate p-values is addressed. Extensive Monte-Carlo sim¬ 
ulations show the adequate finite-sample behavior of the studied tests for a large 
number of realistic data generating scenarios. Illustrations on several environmental 
datasets conclude the work. 
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1 Introduction 


The block maximum method put forward in the seminal monograph of Gumb el (119581 ) 
is frequently applied in environmental sciences to analyse extremes of a given series of 
observations. It consists of dividing the initial high frequency observations into blocks 
of equal size and forming a sample of block maxima. A typical example is the study 
of annual or seasonal maxima of temperatures or precipitation from daily observations. 
The distribution of the block maxima is classically modeled using the generalized extreme 
value (GEV) distribution. Under rather broad conditions on the under lying initial obser¬ 
vations, the extremal types theorem (see, e.g., de Haan and Ferreiral. 20061) states that, 
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as the block size increases to infinity, the only possible non-degenerated limit of affinely 
normalized block maxima is the GEV. The cumulative distribution function (c.d.f.) of 
this three-parameter distribution is 


Gn,<r,zi X ) 








J 


if f 7^ 0, 

if e = o. 


(i.i) 


for x G M, where /i 6 1, ff > 0 and ( 6l are the so-called location, scale and shape 
parameters, respectively, and with the notation (•)+ = max(-,0). 


The adequacy of the GEV distribution as a model for a sample of block maxima 
is predicted by the theory when the underlying high frequency observations are identi¬ 
cally distributed and independent, or satisfy some suitable mixing conditions (see, e.g., 


Leadbetter et ah, 1983 1). and the block size is “sufficiently large”. Sometimes maxima are 


directly observed, making it impossible to question the stationarity of the initial series. 
Furthermore, the block maxima method may be applied even when the initial series is 
clearly not identically distributed. For instance, when the original observations exhibit 
seasonality, it suffices to take the block size to be a multiple of the season length. 


When the block size is large, one of the least questionable assumptions might be 
that of independence among block maxima. Before fitting a GEV distribution to the 
resulting extremes, two questions then arise naturally: (i) Is the series of independent 
maxima identically distributed? (ii) Is the block size large enough to warrant a good fit 
by the GEV distribution? The aim of this work is to focus on the first question without 
assuming that the answer to the second question is positive. Specifically, we consider that 
we have at hand a sample of independent block maxima Xi,, X n and our objective is 
to develop tests for 


Hq : 3 F such that Afi,..., X n have c.d.f. F 


( 1 . 2 ) 


against alternatives involving the non-constancy of the c.d.f. In other words, we aim at 
developing tests for detecting changes in the distribution of independent block maxima. 

Under the additional assumption that each X t is GEV distributed, a natural approach 
for testing H 0 is to apply the maximum likelihood based t heor y dev eloped in the first 
chapter of the well-known monograph of iGsorgo and Horvath! (1199711 on change-point 
analys is. The latter approach was considered for instance in Jaruskova and Rencova 
(2008) but not without significant practical difficulties related to maximum likelihood es- 
ti mation of the GEV parameters from small samples (see also for instance the discussions 


m 


Hosking et ah . 19851: Diebolt et ah . 2008 b 


The approach that we consider is based on yrobability weighted, moments (PWM) (see, 




Greenwood et 


2008; Guillou et al 


2009; 


1979: Hosking et ah, 1985 


Hosk ing and Wallis, 1987; Diebolt et al. 


Ferreira and de Haan, 2015) and has the advantage of giving 


rise to tests that can hold their level asymptotically even if the distributions of the block 
maxima Aj are not GEV. To illustrate that point, in our simulations, some scenarios 
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involved block maxima obtained from blocks of size 1, 5, 10 and 50. The later sen¬ 
tence actually implies that the tests can be applied outside of the block maxima setting, 
although, by construction, this is less natural. These aspects will be clarified in the 
forthcoming sections. 


The rest of the paper is organized as follows. In Section [2j we recall the PWM method 
for estimating the parameters of the GEV due to iHosking et al] (119851 ). propose related 
cumulative sum (CUSUM) statistics for change-point detection, establish the limiting 
null distributions of the latter and address the computation of approximate p-values. The 
third section proposes alternative tests in the framework of the generalized PWM method 
for estimating the GEV due to Diebolt et ah (120081 ). Simulation results are partially 


reported in the fourth section and show the adequate behavior of the tests for many 
realistic data generating scenarios. The last section presents several fully reproducible 
illustrations on environmental datasets. 


All proofs are deferred to a sequence of appendices. Appendices |B] and O hi particular 
establish necessary intermediate results that might be of independent interest, namely 
the weak convergence of the sequential weighted uniform empirical process and of certain 
sequential empirical processes constructed from PWM estimators. The paper is also 
accompanied by supplementary material available upon request from the authors and 
mostly consisting of simulation results that are not reported in Section [4j Finally, note 
that the stud i ed te sts for change-point detection are implemented in the package npcp 
(IKoiadinovid. 1201511 for the R statistical environment ()R Gore Team! . 120161) . 


2 The PWM method and related statistics 


The statistics for testing H 0 in (11.21) th at we propose ar e related to the method of PWM for 
estimating the GEV parameters (IHosking et all . 1985)- After briefly recalling the method 
and suggesting related statistics for change-point detection, we study the limiting null 
distributions of the latter and address the practical computation of approximate p-values. 


2.1 Estimation of the GEV using PWM 


A general definition of PWM can be found for instance in Hosking et al. (119851 Section 2) 
(see also Appendix ED- Let X be a random variable with c.d.f. F and consider the func¬ 
tions vi(x) = 1, zz 2 (x) = x and v 3 {x) = a; 2 , x G [0,1]. The quantities j3i = E[Xvi{F(X)}], 
i e {1, 2, 3}, are the three (probability weighted) moments used in the PWM method for 
estimating the GEV parameters. 

Assume that H 0 in (11.21) holds with F the c.d.f. of the GEV given in (II.ip . The 
method exploits the following relationship between j3i, j3 2 and f3 3 and the parameters /i, 
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cr and £ of the GEV. For a > 0, £ < 1 and £ 7 ^ 0, there holds 


f A = M-|{i-r(i-£)}, 

20 2 - 0i = |t(1-£)(2*-1), 

3/3 3 - A = 3^ — 1 

l 2 «-l’ 


( 2 . 1 ) 


where T is the Gamma function. Using continuity arguments, some thought reveals that 
the constraint £ 7 ^ 0 is actually unnecessary. Estimators of /x, a and £ are then the 
solution of the above system of equations in which 0 i, 0 2 and 0 3 are replaced by their 
estimators. 


Two ways to estimate 0i, 0 2 and 0 3 were considered by I Hosk in g et al, ( 1985 ). The 
first approach consists of first estimating F by a slight modification of the empirical c.d.f. 
of the sample, and then estimating the PWM by suitable arithmetic means. For reasons 
that shall become clear later in this section, we present the resulting estimators when 
based on a subsample Xk, ■ ■ ■ ,Xi, 1 < k < l < n, of the available data. The c.d.f. F is 
estimated by 


F k: ,(x) = 


1 


l — k + 1 


^ 1 (Xj < x) + 7 > , i 6 l, 


( 2 . 2 ) 


. j=k 


with the convention that TV = 0 for k > l and where 7 is a constant whose value 
suggested by [Hosk ing et al.l (119858) is -0.35. A sensible estimator of 0* is then 


Pi,k:l 


l — k + 1 ^ 

j=k 


Y, XMF tJ (X,)} 


(2.3) 


with the convention that 0 *^.; = 0 for k > l. 

A second equally natural possibility to estimate 0i, 0 2 and 0 3 is to con sider the 
asymptotically equivalent estimators & 2 ,i:n, & 3 ,i :n proposed by Landwehr et al. (119791 ) 
and defined, for a subsample Xk, ..., Xi, 1 < k < l < n, by 


bi,k:l 


1 y^ rKi{(j-fc + l)-™} 

l ~ k + l ht -k + l)-m} 


X, 


(j),k:h 


(2.4) 


where X( k ),k-.i < • ■ ■ < X^^a are the order statistics computed from X k , ■ ■ ■, X h and with 
the convention that = 0 if l — k < i. The latter estimators are unbiased (they are 
related to [/-statistics). 

, Section 4) claimed that the 0^,/ are slightly superior to the 
in terms of finite-sample performance. While this seems true when X \,..., X n are i.i.d. 
from a GEV with location zero and scale one, it does not seem to necessarily hold for 
arbitrary location and scale parameters. 

Let g /t , g a and be the components of the map from R 3 to R 3 implicitly defined 
by (12.ip that transforms f3 = (01,02,03) into (/q cr, £) (that is, /j = g^f 3), a = g a {(3) and 


H oskin g et ahj ( 1985 
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£ = g^(/3)), and let us first focus on the function without assuming anymore that F is 
the c.d.f. of the GEV but only that F is continuous and that Aj A and A exist. Under 
the constraint £ < 1 , it can be easily verified that the solution g^(/3) of the third equation 
in ( 12 .ip exists and is unique if and only if 1 < (3A — /5i)/(2/3 2 — A) < 2 . The latter is 
guaranteed by the following result proved in Appendix [~Y1 

Proposition 2.1. Assume that F is continuous and that A? A an d A exist. Then, 

2A - A > 0 , 3A - 2A > 0 and -A + 4A - 3A > 0 . 


Given the equivalence 

r 2A-A>0, f 2A - A > o, 

{ 3A-2A>0, « l 3A-A>0, 

[ —A + 4A — 3A > 0 , [ l < (3A — A)/(2A — A) < 2 , 

it is natural to consider g% as a function defined on the convex subset 

T>£ = {x e M 3 : 2x2 — oq > 0, 3x 3 — 2x 2 > 0, —x\ + 4x 2 — Sx 3 > 0}. (2.5) 


In addition, it is easy to see that, since (3 e necessarily g%{(3) < 1 and g a (f3) > 0, 
which is fully in accordance with the constraints under which the system ( 12 . lj) was derived. 
Using the fact that x e->- x/[2 x — 1) can be extended by continuity at 0 and is strictly 
positive, it is also natural to consider g a and g lt as functions defined on 


When Aj A and A are estim ated u sing the unbiased estimators 6 1 > 1 :n ,A,i:reA 3 ,i:n 
dehned analogously to (12.41) . [Hosking et ai- (119851 . Appendix B) showed that A :n = 
A,i:n) A,i:n) £ 2 ^ f° r all 71 — 3, where T>£ is dehned in ( 12 .5p . which implies that 
g^ipi-n) < 1 and g a (bi :n ) > 0 for all n > 3. The latter desirable property is referred to 
as the feasibility criterion of the PWM method. When 0i :n = (Ap :n , A,i :n , A,i:n) is used 
instead of bi :n , an analogue result is not available. However, provided E(|Xl|) < oo, a 
consequence of the strong law of large numbers and the Glivenko-Cantelli lemma is that 
(3i :n converges almost surely to (3 (see Appendix [D]) . The latter immediately implies that 

P(/An ^ — > 1 . 


2.2 Test statistics 


The above developments suggest to study tests for change-point detection based on the 
three CUSUM like statistics 


Snn = 


k(n — k ). 


'g,n - , max , - 370 —!(/3l:fc e l:n £ A/A:fc) ~ AA+l:n) , 9 £ {g„, 9a, 9$}, 

1 </e<n—1 TL ' 


( 2 . 6 ) 


or on the three alternative ones 


fc(n — A;) 
A,n = max -^— 

3<k<n-3 n 6 ' 2 


9(bl:k) 9(pk+l:n) 


9 ^ {S'm’ .9?}- 


(2.7) 


In (USD, /3nfc = (A,i:fc, A,i:fe, A,life) and /A-+i :n = (A,fc+i:ra, A,fc+i:n^A,fc+i:n) are dehned 
according to (I2.3j) . and V^ is given in (j2.5j) . while in (I2.7j) . 61 ,* = (A,i:fc,Ai:fcA 3 ,i:fc) and 
bk+i-.n = (bi,k+i:n,i> 2 ,k+i:n,h,k+i:n) are dehned according to Q- 
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For instance, under Hq in (ll.2[) with F the c.d.f. of the GEV, S g ^ n is the maximum over 
k G {1,. .., n — 1} of the normalized absolute difference between g^{(3i±), the estimator 
of the shape parameter £ computed from Ah,, X k , and g^($ k +i :n ), the same estimator 
computed from X k+ i,... , X n . The role of the indicator l(/3 1:fc G T>£,/3 k+ i :n G T>^) is to 
allow the previous evaluations of g^ only if (3 1:k G V^ and 0 k +i:n G The coefficient 
k(n — k)/n 3//2 in front of each absolute difference in (12.611 is the classical normalizing 
term in the CUSUM approach ensuring that, under additional suitable assumptions, 
S gen converges in distribution under the null. The interpretation of T„ n is similar. The 


reason for the disappearance of the indicators in i ts exp ression is that, as discussed in 
Section [2TTT b\. k G T>^ for all k > 3 ( Hosking et ah 


1 985 . Appendix B). 


If each block maxima in the sample X\,... ,X n is GEV distributed, S g( _ n and T< 


95 >' 


(resp. S g<j) n and T g<j n , S g n and T g ^, n ) are thus test statistics particularly sensitive to 
changes in the shape parameter £ (resp. the scale parameter cr, the location parameter /i) 
of the GEV. It is however important to keep in mind that the proposed tests are actually 
nonparametric in the sense that they remain potentially meaningful even if the block 
maxima are not GEV distributed. The latter follows from the fact A.i:n = Ann is simply 
the sample mean, while, for instance from Hosking et ahl ( 1985 . page 252), we have that 
2 & 2 ,i:n — &i,i:n is actually Giiii’s mean difference and, under H 0 in ( 1 1 . 2 j) . is thus an estimator 
of scale as measured by 2A _ A = E(|Ai —A 2 |)/2. Consequently, the tests based on T g ^^ n 
and S g n (resp. T ga , n and S g<T _ n ) can in general be interpreted as tests for change-point 
detection particularly sensitive to changes in the location (resp. scale) of the distribution. 
Similarly, the tests based on T g ^ n and S^.n can be interpreted in general as tests for 
change-point detection particularly sensitive to changes in the upper tail of F. Indeed, 
the ratio (3 A — A)/ (2 A — A) appearing in the expression of g^ can be rewritten as 


E{max(A 1 , X 2 , X 3 ) - (Ah + X 2 + X 3 )/3} 

TVmnYl (Y, - (Y, 4- 


and is thus clearly invariant to the replacement of the A \ by cX t + d if c > 0. Some 
thought reveals that this ratio captures some information on the upper tail of F. 

The statement that the tests are actually nonparametric will also become clearer 
in view of the theoretical results established in the next two subsections and will be 
illustrated in the Monte Carlo experiments partially reported in Section 


2.3 Limiting null distributions of S g , n and T g . n 

We shall now state conditions under which S g . n in ( 12 . 61 ) and T gn in ( 12 . 71 ) converge in 
distribution under the null. 

Condition 2.2. The random variables Xi,... ,X n are i.i.d. from a continuous distri¬ 
bution with c.d.f. F such that E(A^) < oo and there exists a G [0,1/2) such that 
sup xeR \H a (x)\ < oo, where H a (x) = x[F(x){ 1 — F(x)}] a , iGl. 

Condition 12.21 is for instance satished when F is the c.d.f. of the GEV in ( 11 . 11 ) 
with £ < 1/2, that is, when it has finite variance. In addition, it can be verified that 
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when both max(Xi,0) and — min(Ab,0) have an ultimately monotone density and are 
in the maximum domain of attraction of the GEV with shape parameter £ < 1/2, then 
su Pa;eR < OO for any a G (£,1/2) (see Section 3 in the supplementary material). 

In the sequel, the first-order partial derivatives of g G {.9/j, ,9<x> fk\ be denoted 
by dig, d 2 g and d 3 g, respectively. Recall furthermore that n\{x) = 1, u 2 (x) = x and 
v 3 (x) = x 2 , x G [0,1]. 

The following result is proved in Appendix iDl 

Proposition 2.3 (Limiting null distributions of Sg, n and T gyTl ). Under Condition \2.2l for 
any g G {g p , g a , g^}, Sg^ n —T gn = op(l) and S 9jn converges in distribution to a g sup sg [ 01 ] |HJ(s) — 
sU(l)|, where U is a standard Brownian motion on [0,1] and 

3 

v g = Y d i9(P)d3g(P)Cov(Y Vi ,Y v .), (2.8) 

i,j=l 

with, for any v G {zq, v 2 , u 3 }, 

Y u = Xii/{F(Xi)} + [ xu'{F(x)}l(Xi < x)dF(x). (2.9) 

Jr 

2.4 Computation of approximate p-values 

Proposition 12.31 suggests to base the computation of approximate p-values for S g>n and 
T g n , g G {g p , g a , g^}, on their estimated asymptotic null distributions. Indeed, the dis¬ 
tribution of the supremum of a Brownian bridge, known as the Kolmogorov distribution, 
can be approximated very well in practice. It therefore simply remains to estimate the 
unknown variance a 2 in (12.81) . 

To do so, it is necessary to estimate the first-order partial derivatives of g at /3 and the 
covariances Cov(Y/., Y Vj ), i,j G {1, 2, 3}, where Y v . is dehned in (12.91) . Natural estimators 
of the former are <9ig(/3i :n ), . .., d p g(/3i :n ) or dig{hi, n ), ..., d p g{bi, n ), while estimators of 
the latter can be based on the pseudo-observations 

1 n 

Y Vli , n = XiViFmiXi)} + -Y X ^'{ F l:n{X ] )}l{X l < Xj), 

U 3=1 

for i G {1,... ,n} and v G {zq, z/ 2 , v 3 }, where F x . n is dehned in (I2.2j) . Indeed, a sensible 
estimator Co v n (Y Vi ,Y Uj ) of Cov(Y Vi ,Y Vj ), i,j G {1,2,3}, under H 0 is then simply the 
sample covariance computed from (Y Ui)hn ,Y Uj ^ n ),... ,(Y ViiUtn , Y^. ,„,„). 

The proof of the following result is omitted for the sake of brevity. 

Proposition 2.4. Under Condition \2. 21 for any g G {g g , g a , g^}, 

p 

Vg,n = Y d i9(Pi:n)d j g((3i :n )Cov n (Y „ i , Y v .) (2.10) 

i,j =1 
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and 


& ln = Y1 ^(^l:n)^(bl:n)C0V n (y i/ .,y i/ . 

i,j =1 


converge almost surely to a 2 in 


As we continue, we shall thus compute approximate p-values for Sg t71 as 1 —F^S^/cr^n) 
and approximate p-values for T g;n as 1 — FK(T g ^ n /a gjn ), where Fk is the c.d.f. of the Kol¬ 
mogorov distribution. As classically done, we approached Fk by the c.d.f. of the statistic 
of the Kolmogorov-Smirnov goodness-of-ht test for a simple hypothesis. From a practi¬ 
cal perspective, we used the function pkolmogorovlx given in the code of the R function 
ks.test. 


An even more important practical issue is the computation of S g n and T g n for g equal 
to <? M , g a and g^. The numerical solving of the system (12.11) can actually be avoided by 
using the following accurate approximations proposed in 


lHosking~ 


et al. (1985): 


£ = 9 e(P) = f ° fdP) with 


& = 9 a{P) 


(2(3 V2 -p vi )j 

r(i-0(2f-i) 


h = 9n(P) = - r(i 


f(z) 


—7.8590a; - 2.9554a; 2 , 


M x i,x 2 ,x 3 ) 


2 x 2 — x\ log 2 
3aq — x\ log 3’ 


1 


-«}• 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


We also used the above approximations in combination with symbolic differentiation to 
obtain approximations of the first-order partial derivatives of g fl , g a and g^ necessary for 
computing a 2 n and d g n in Proposition 12.41 


3 The GPWM method and related statistics 


Alternative test statistics could yet be based on the generalization of the PWM method 
for estimating the parameters of the GEV proposed by Diebolt et ah (2008)- Under 
the assumption that H 0 in (II. 2\i holds with F the c.d.f. of the GEV distribution, the 
estimation method involves the quantities = E[Xu i {F(X)}], i G {1, 2, 3}, with Ui(x) = 
—a: log ay u 2 (x) = x(logx ) 2 and v?.[x) = — x 2 log x, x G [0,1]. The moments /?i, (3 2 and f3 3 
are called generalized PWM in Di ebolt et al.i 020 0 8 ) as t hey are not PWM in the classical 
sense considered for instance in Greenwood et ah ( 19791 ) because of the definitions of the 
functions zq, v 2 and zq. 

The relationship between the GEV parameters /a, a and £ and the GPWM f3i, (3 2 and 
/3 3 is given by the following system of equations, which holds for a > 0, £ < 2 and £ 7 ^ 0: 


' 4/3 1 = M -|{l-2T(2-£)}, 

< A-/3 2 = ^r( 2 -£), 

2Q01-&) . e 

, A-9/4ft l-(3/2)V 

8 


(3.1) 




































Again, the constraint £ 7 ^ 0 can be removed using continuity arguments. Estimators of fi, 
a and £ are then the solution of the above system in which f3i, fa and fa are replaced by 
suitable estimators. A sensible choice is to consider fa t i :n , $ 2 ,i :n and fa t i :n , respectively, 
defined analogously to ( 12 . 3p . 

Let hh a and h ^ be the components of the map from IR 3 to M 3 implicitly defined 
by (13.11) that transforms (3 = (/3i, /3 2 , /3 3 ) into (/ 1 , cr, £). In practice, to evaluate whether 
/3i:n = fa.i-.n, fai-.n) £ IR 3 belongs to the domain of definition D h of h £, h a and h 

we shall attempt to solve numerically (13.11) and will conclude that fa, n G T>y , if and only 
if fa(/3 1:n ) < 2, h a ((3 i :ri ) > 0 and h 9 (f3 1:n ) G IR. 

The above ingredients lead to the following generalization of the statistics Sg }Tl in (12.61) : 


S h<n = max 

l<fc<n—1 


k(n — k ). 


n 


3/2 


l(/3l :k £ '^hi /3fc+l:n ^h) i h G { hg , h a , } . 

(3.2) 


Should the block maxima A 7 "* be GEV distributed, Sh ( , n (resp. Sh a , n , Sh^^n) is a test 
statistic particularly sensitive to changes in the shape (resp. scale, location) parameter. 


We conjecture that the asymptotics of the statistics in (13. 2j) under H 0 in (II . 2 j) take 
the same form as those of the statistics S 9tn considered in the previous section, the only 
difference being the dehnitions of the functions iq, u 2 and v 3 . Let a\ be the analogue 
of (12.81) . Specifically, we conjecture that, under the null, Sh, n converges in distribution to 
ah sup sG [ 01 ] |U(s) — sU(l)|, where U is a standard Brownian motion on [0,1]. By means 
of a sequential extension of the method of proof used in Diebolt et ah ( 20081) . it might be 
possible to show the aforementioned weak convergence when F is the c.d.f. of the GEV 
and h is either fa, h a or h 9 . 


To carry out the tests, consider the analogue a\ n of ( 12 . 10 p . Approximate p-values 
for Sh.n can be computed as 1 — FK(Sh, n /&h,n ), where Fk is the c.d.f. of the Kolmogorov 
distribution. 


From a practical perspective, as in the previous section, the numerical solving of 
system (13.11) can be avoided. Starting from the fact that the function x (-)■ {x/(l — 
(3/2) a: )} log ) 3//2 ) is almost linear on [—1,1], we propose the following accurate approxima¬ 
tions of fa, fa and /i /t : 


f = = f 0 with 


a = fa(f3) 


2 3 -^(A-/g 2 ) 

r(2 -0 


f(x) 


1.442853 - (- x )°- 4 °5465i 
0.1183375 


h( x 1,X2,X 3 ) 

K 


2 (xi - x 2 ) 

Q ) 

Xi - fx 3 


n = M/3) = 4A + |{i - 2 s r (2 - {)}. 


(3.3) 


(3.4) 

(3.5) 


As in Section 12.41 approximations of the partial derivatives of he, h a and h fJ which are 
necessary for computing cf 2 n , are obtained by differentiating , fa and h lt . 
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4 Monte-Carlo experiments 


Monte-Carlo experiments were carried out in order to investigate the finite-sample be¬ 
havior of the tests studied in the previous sections. Nine statistics were considered: S g n 
in (12.61) and T gn in (12.71) with g being respectively in (12.lip . g a in (12. 12j) and g g 
in (12.131) . and the statistic Sh, n in (13.21) with h being respectively h^ in (13.3j) . h a in (13.4j) 
and hfj, in ([375]) . Recall that, under the conditions of Proposition 12.31 the statistics S g>n 
and T g>n are asymptotically equivalent under the null, and that the functions g%, g a , g g 
(resp. fi£, h a , hg) are related to the PWM (resp. GPWM) method for estimating the GEV 
parameters summarized in Section I27T1 (resp. Section [3]). When computing the statistics, 
the following asymptotically negligible change was made: given a small integer r > 1 , 
maxima in (12.61) . (12.71) and (13.2[) were taken over the set {r,... ,n — r} instead of the set 
{1,..., n — 1} or {3,..., n — 3}. The role of the integer r is to exclude subsamples of 
cardinality strictly smaller than r. In our simulations, we used the value r = 10 and, 
in all the considered settings, encountered hardly any sample for which the indicators 
in (12.61) and (|3. 2 j) turned to be zero for some value of k. 

We started our experiments by investigating the empirical levels of the nine tests for 
samples of size n G {50,100,200,400} generated from the GEV(/qcr,£). To estimate the 
power of the tests, 1000 samples were generated under each data generating scenario and 
all the tests were carried out at the 5% significance level. Our main finding was that all 
tests had a tendency of being too liberal when \g/a\ was large. This prompted us to run 
the tests based on S gn in (12.61) (resp. T g>n in (12 . 71) . S h n in (13.21) 1 on the data translated by 
-<7 m (/ 3i:n) (resp. -g g {b 1:n ), -/i M (/3i :n )). Since g g ((3 1:n ) ^ <7 M (/3) and g^(b 1:n ) -^4 g M (/3) 
under Condition 12.21 the asymptotic results stated in Section 12.31 remain valid when the 
tests based on the statistics S g _ n and T gn are applied on the translated data. 

The adopted translation of the data led to a significant improvement of the empirical 
levels. Additional improvement of the levels of the test based on T gatn (resp. Tg (tU ) was 
obtained by multiplying cr^ n (resp. <rj n ) in ( 12 . 101 ) by {n + 10 )/n (resp. (n + 20 )/n). 
A last asymptotically negligible change was carried out for the tests based on Sh, n - the 
statistics given in (13.2j) were computed with 7 in (12.2j) equal to zero. 

As the tests based on S g _ n in (12.61) turned out to be consistently worse-behaved than 
those based on T g n in (12.71) . for the sake of brevity, we do not report any results for the 
former in the forthcoming tables. 

Table [T] gives the empirical levels of the six remaining tests for samples generated 
from the GEV(/i,l,£) for g = 0 and £ G { —1, —0.9,..., 1,1.1}. Analogue tables for 
g G { — 100, —10,10,100} can be found in Section 1.1 of the supplementary material (see 
Tables 1-4 therein). For the sample sizes under consideration, the tests based on T g ^ n 
appear to hold their level well for — 1 < £ < 0.4, that is, as long as the variance of the 
GEV is finite. The tests based on 5} )Tl are overall too conservative for — 1 < £ < 0.3. 
Their empirical levels are closer to the 5% nominal level for 0.4 < £ < 0.8 and they start 
being too liberal for £ > 0.9. 

A similar experiment was carried out for samples of independent block maxima ob¬ 
tained from blocks of size b G {1,5,10,50} from the generalized Pareto distribution 
(GPD) with location parameter equal to 0, scale parameter equal to 1 and shape pa- 
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rameter equal to £ (abbreviated as GPD(0,1,£) as we continue). The experiment was 
designed having in mind the fact that the GPD with shape parameter £ is in the max¬ 
imum domain of attraction of the GEV with shape parameter £. The results for 6=1 
are reported in Table [21 The results for 6 G {5,10, 50} are reported in Section 1.2 of the 
supplementary material (see Tables 5-7 therein). Another variant of this experiment was 
carried out by using the absolute value of the standard Student t distribution with l/£ 
degrees of freedom instead of the GPD(0,1,£). Indeed, the absolute value of the standard 
Student t distribution with l/£ degrees of freedom, £ > 0, is in the maximum domain of 
attraction of the GEV with shape parameter £. The results are reported in Section 1.3 
of the supplementary material (see Tables 8-11 therein). The conclusions in terms of £ 
are overall the same as for the first experiment: the tests based on T gn behave mostly 
adequately when £ < 0.4 and become too liberal for £ > 0.5; the tests based on Sh, n are 
overall too conservative for £ < 0.3 and are too liberal for £ > 0.9. 


[Table 1 about here.] 
[Table 2 about here.] 


Additional experiments under the null (whose results are reported in Table 12 of the 
supplementary material) were carried out to assess the empirical levels of the tests for 
samples generated from various normal and exponential distributions. No test turned out 
to be too liberal. 


The next experiments focused on the power of the tests. To ease reading of the 
forthcoming tables, among the six tests proposed in the paper, the rejection percentages 
of those that are expected to be sensitive to the alternative under consideration are 
colored in light gray. 


We started by assessing the power of the tests in the case of a change in the shape 
parameter of the GEV. First, for t G {0.25,0.5,0.75}, samples of size n G {100,200} 
were generated such that the \_nt\ first observations were from a GEV(/ql,-0.4) and 
the n — [nt\ last observations were from a GEV(/r,l,£), for £ G {—0.2, 0, 0.2, 0.4} and 
/i G {—100,0,100}. The results for /i = 0 are reported in Table El while those for 
/i G {—100,100} are given in Tables 13 and 14 of the supplementary material. The 
column F n gives the rejection percentages obtained using the nonparametric rank- based 
test for change-point detection based on empirical c.d.f.s studied in Holmes et al. (2Q13|) 
and implemented in the function cpTestFn of the R package npcp. The columns x n and 
report the empirical powers of CUSUM tests designed to be particularly sensitive to 
changes in the expectat ion (see [ Phillips!. Il987l. and the references therein) and the vari¬ 


ance (see, e.g. 


Bucher and Koiadinovic . 201 til ), respectively. The former (resp. latter) 


is 


implemented in the function cpTestMean (resp. cpTestU) of the R package npcp. Ta¬ 
bles [4] and [5] report similar results but for higher £ values. Overall, as could have been 
expected, the tests based on T g ^ n , T ga , and n have hardly any power against 

such alternatives. Roughly speaking, the test based on T^ >ri is more (resp. less) pow¬ 
erful than the one based on S^ n when the largest shape parameter is smaller or equal 
than 0.2 (resp. greater or equal than 0.4). Both tes ts are ove rall mo re p owerful than 
the general-purpose nonparametric test considered in [Holmes e t al. (120131) and the two 
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CUSUM tests designed to be particularly sensitive to changes in the expectation and the 
variance, respectively. 


[Table 3 about here.] 
[Table 4 about here.] 
[Table 5 about here.] 


Similar experiments were used to assess the power of the tests when the scale (resp. 
location) parameter increases from 0.5 to 1 (resp. from 0 to 0.5). The corresponding 
rejection percentages are reported in Tables [6] and [7] respectively (additional results are 
available in Sections 2.2 and 2.3 of the supplementary material). For changes in the 
scale parameter, we see that the test based on Tg a>n is more (resp. less) powerful than 
the test based on for £ G {—0.4,0} (resp. £ G {0.4, 0.8}), the two tests becoming 
more powerful than the test “s^” for £ > 0. Again, as expected, the remaining PWM 
or GPWM tests have little power against such alternatives. For changes in the location 
parameter, the test based on Tg n is more (resp. less) powerful than the test based on 
for £ G {—0.4,0} (resp. £ = 0.8). This time the rank-based test of lHolmes et ah 


Sl 


(120131 ) is more powerful than the two aforementioned tests, although the difference is 
small for £ G {0, 0.4, 0.8}. 


[Table 6 about here.] 

[Table 7 about here.] 

In a last experiment, we investigated the power of the tests for samples generated from 
normal distributions with different variances or expectations. The results are reported 
in Tables 23 and 24 of the supplementary material. In the case of a change in variance 
(resp. expectation), the proposed tests designed to be sensitive to a change in scale (resp. 
location) appear overall slightly less powerful than the (resp. “x n ”) tests. 


5 Illustrations 

The results of the finite-sample experiments partially reported in the previous section 
suggest that the overall behavior of the tests based on the statistics T fhn in (12.71) is better 
than that of the tests based on the statistics Sh, n hr (13.2[) when the distributions of the 
block maxima have finite variances. When they have an infinite variance but a finite 
expectation, the tests based on the statistics Sh, n should however clearly be preferred as 
those based on T g , n do not hold their level anymore. 

When dealing with climate data, the assumption of finite variances is reasonable for 
temperatures, precipitation or winds, which is why we focus on the tests based on the 
statistics T g>n in the forthcoming applications. 
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To illustrate the pratical use of the three selected tests, we ap plied them on some of 
the environmental datasets available in the R packages evd (IStenhensonl . 1200211 . extRemes 
( Gilleland and Katz . 2011 1. ismev ( Heffernan and Stephenson. . 2014]) and ClusterMax 


(Bernard et all 2 0 13). The observations consist of block maxima of wind speeds, tem¬ 


peratures, precipitation or sea levels, and can reasonably be considered to be independent. 
The interest reader can easily plot the data after installing the aforementioned R packages. 


The main practical difficulty that prevented a direct application of the tests is that all 
the datasets contain ties due to a limited precision of the underlying measurement instru¬ 
ments. To address this issue, we proceeded as follows. For every dataset, we computed 
the smallest absolute difference d between any two distinct observations. Assuming that 
d is a good estimate of the precision of the measurements, we added a random component 
uniformly distributed over (0, d) to every observation. For every dataset, 1000 different 
de-tied samples were generated with the hope that at least one of these is close to the 
true unobserved sample, and the tests were applied to each such de-tied sample. The 
minimum and the maximum of the obtained p-values are reported in Table [HJ The mini¬ 
mum and the maximum of the estimates <7 M (&i/,i:ri), i:n ) and g^{b Vt i :n ) are provided 

as well. Under H 0 in (II. 2p with F the c.d.f. of the GEV, the three latter quantities are 
estimates of the corresponding location, scale and shape parameters obtained using the 
PWM method summarized in Section l2Tl 


[Table 8 about here.] 

Small maximal p-values provide evidence against H 0 in (11.21) . Under the additional 
assumption of GEV distributed block maxima, the latter can be interpreted as evidence 
of change in the location, scale or shape parameter. Ignoring a necessary adjustment 
of the p-values due to multiple testing, the following conclusions could for instance be 
drawn from Table [SJ 

• there is very strong evidence of a change in the location parameter for the annual 
minimal temperature in the HEAT dataset, and strong evidence of the same nature 
for the annual maximal temperature in the HEAT dataset and for the annual maximal 
sea levels in the fremantle dataset, 

• there is weak (resp. some) evidence of change in the shape parameter for the weekly 
maxima of hourly precipitation measured at Saint-Lizier (resp. Nevers), France, 
during the Fall season between 1993 and 2011, 

• there is strong evidence of a change in the scale parameter and some evidence of 
a change in the location parameter for the weekly maxima of hourly precipitation 
measured at Niort, France, during the Fall season between 1993 and 2011. 

Notice that the tests, as implemented in R package npep, are fast. Obtaining the 
results reported in Table [8] took less than 30 seconds on one 2.7 GHz processor. 

Let us end this section with a practical remark. To test Hq in (11.21) from a single sample 
of block maxima, three tests, based on the statistics Tg >n in (12.71) or on the statistics Sh, n 
in (13.21) . are to be used. As noted by a referee, such an approach, although interesting 
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from a data analysis perspective, lacks a simple decision rule regarding the rejection or 
not of Hq. The later could be the subject of future research. A simple yet conservative 
approach would be to use a Bonferroni correction and reject Hq as soon one p-value is 
below a third of the chosen significance level. 
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A Proof of Proposition 12.1 

Proof. Some thought reveals that 

A = E[A{F(A)}* _1 ] = r 1 E{max(A 1 ,..., A»)}, i e {1, 2, 3}, 

where Ad,..., Ad are independent copies of X. 

The first inequality is obtained by noticing that max(A 1 , X 2 ) > (Ad + Xf)/2 almost 
surely since P(Ad = X 2 ) = 0 by continuity of F, and by taking expectations on both sides. 
The second inequality follows from the fact that 3 max(Ad, Ad, X 3 ) > max(Xi, Ad) + 
max(Ad, AG) + max(X 2 , X 3 ) almost surely. The last inequality is a consequence of the 
fact that 

2 max(Ad, Ad) + 2 max(Ad, X 3 ) + 2 max(X 2 , X 3 ) > 3 max(Ad, X 2 , X 3 ) + Ad + Ad + Ad 
almost surely. ■ 


B Weak convergence of the sequential weighted uni¬ 
form empirical process 


The purpose of this appendix is to derive the weak convergence of a sequential version of 
the weighted uniform empirical process. The latter result, that may be of independent 
interest, is used in the proof of Proposition lC.21 itself necessary for the proof, given in Ap¬ 
pendix [D] of Proposition 12.31 The forthcoming developments rely heavily on Appendix G 
of Genest and Segersl (2009). 


Let U\,... ,U n be an i.i.d. sample from the standard uniform distribution. For any 
integers 1 < k < l < n, let Gk-.i be the empirical c.d.f. computed from the sample 
Uk,... ,Ui with the convention that Gk-i = 0, for all k > l. Then, define the two-sided 
sequential uniform empirical process by 


[nt\ 


(s, t, u) = — {l(Ui < u) -u} = v / ^A n (s,t){G Lns j + i :Ln q(u) -u}, (B.l) 


z=[nsj+l 
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for (s,t,u) G A x [0,1], where A = {(s, t) G [0, l] 2 : s < t} and \ n (s,t) = ([nt\ — 
|_nsj)/n, (s,t) £ A. The process B n (( h •) is the standard sequential u niform empirical 
process and it is well-known (see, e.g. Jvan der Vaart and Wellner . 2000 . Theorem 2.12.1) 
that B„(0, •, •) B(0, •, •) in £°°([0, l] 2 ), where B(0,-,-) is a centered Gaussian process 

with continuous trajectories, sometimes called a Kieffer-Miiller process , whose covariance 
function is given by 

Cov{B(0, s, u), B(0, t,v)} = (s A t)(u A v — uv ), (s, t, u, v) £ [0, l] 4 . 

Since, for any (s,t) £ A, M n (s,t, •) = B n (0, t, •) — B n (0,s, •), we immediately obtain 
from the continuous mapping that B n -w B in £°°(A x [0,1]), where, for any (s,t) £ A, 
B(s, t, •) = B(0, t, •) — B(0, s, •). 


Analogously to Genest and Segers (2009), let a > 0 and consider a weighted version 
of the processes B n defined as 




W a , n (s,t,u) = { {«(l- w )}*’ lf ^ £ (0, f), 
0, if u £ {0,1}. 


(B.2) 


The following result is then a univariate sequential extension of Theorem G.l of Ge ne st and Segers 

fteoooh . 

Proposition B.l. For any a £ [0,1/2), W Qin W a in £°°(A x [0,1]), where W a is 
a centered Gaussian process with continuous trajectories such that W Q (0, -,u) = 0 for 
u £ {0,1}, 

(s A t) [u A v — uv) 


Cov{W a (0,s,u),W a (0,f,v)} = 

|u(l — w)u(l — v)\ a 

and W a (s, t, •) = W a (0, t, •) — W a (0, s, •) for all (s, t) £ A. 


(s,t,u,v) £ [0,1] x (0,1) , 


Proof. Fix a £ [0,1/2). Analogously to Genest and Segeri ( 2009 . Appendix G), for any 
u £ (0,1), let f u : (0,1) —» M be defined as f u (x) = {n(l — u)}~ a {l(x < u) — u}, 
x £ (0,1), and let T = {f u : u £ (0,1)} U {0}, where 0 is the function vanishing 
everywhere on (0,JJ. Furthermore, let_P be the uniform probability distribution on (0,1). 
Lemma G.2 of Genest and Segers ( 20091) then states that the collection P is P-Donsker 
(see, e.g., van der Vaart and Wellner . 2000 : Kosorok . 2008h . which is equivalent to saying 
that there exists a P-Brownian bridge G such that G n G in £ 00 (J r ), where, for any 
/ G P, G„(/) = Vn(P n / : - Pf ) with 


r i n 

Pf= f(x) dr, V«f=-Tf(U i ), 

A.i] n i,i 


and the convention that Po/ = 0. For any s £ [0,1] and / G P, let Z n (s, /) = 
i/nA n (0, s)(P[ ns j/ — Pf) = y/n\ n (0, s)P[ ns j/, where the second equality follows from 
the fact that P f = 0 for all f £ P. An immediate consequence of Theorem 2.12.1 of 
van der Vaart and Wellner (2000) is then that Z n Z in £°°([0,1] x P), where Z is a 


tight centered Gaussian process with covariance function 

Cov{Z(s, /), Z(t, g)} = (s A t)(Pfg - PfPg) = (s A t)Pfg, s, t £ [0,1], /, g £ P. 

(B.3) 
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To conclude the proof, we proceed analogously to Genest and Segers 1 (2009, Appendix G). 
Let 4> : [0,1] —> T be defined as 


4>{u) = 


f u , if u e (0,1), 

0, if mg {0,1}. 


Then, let T : ^°°([0, 1] x J) G £°°([0, l] 2 ) be defined as T(z)(s,u) = z(s,4>(u )) for all 
0 G f?°°([0, 1] x J -) and (s,u) G [0, l] 2 , and notice that T(Z n ) = W Qin (0, •, •). Indeed, for 
any s G [0,1], if u G (0, 1), 

T(Z n )(s,u) = Z n (s,f u ) = VnA n (0,s)P Lns j/„ = W Q , n (0,s,w), 
and, if u G {0,1}, 

T(Z n )(s,u) = Z n (s,0) = v / ^A n (0,s)Pi„ s j0 = W Q>n (0, s, u). 


The map T being linear and bounded, it is continuous. Hence, by the continuous mapping 
theorem, T(Z n ) T(Z) in £°°([0, l] 2 ). Let W Q (0, •, •) = T(Z). Given the dehnition of T, 
since Z is centered, so is W a (0, •, •), and since the hnite-dimensional distributions of Z are 
Gaussian, so are those of W a (0, ■, •). Furthermore, the covariance function of W a (0, •, •) 
follows immediately from (IB. 3D . It remains to ver ify that the trajectories of W a (0, •, •) 
are continuous. From van der Vaart and Wellner (2000, Chapter 2.12), we know that 
the trajectories of Z are almost surely uniformly continuous with respect to the natural 
semimetric |s — t\ + p(f,g ) on [0,1] x J where 

p 2 (f,g) = E[{G(f)-G(g)} 2 }, f,gef. 


From Genest and Segersl (2009, page 3020), we have that the map </> is continuous with 
respect to the Euclidean metric on [0,1] and the standard deviation se mi m etric p on T. 
As a consequence, the trajectories of (s,u) i— > W Q (0 ,s,u) = Z(s, 0(w)) are almost surely 
continuous on [0, l] 2 . 

Finally, since W a , n (s,t, •) = ~W a>n (Q,t, •) — W ain (0, s, •) for all (s,t) G A, we immedi¬ 
ately obtain from the continuous mapping theorem that W Qjn W Q in f°°(A x [0,1]), 
where, for any (s, t) G A, W a (s, f, •) = W a (0, t, ■) — W a (0, s, •). ■ 


C Weak convergence of sequential PWM processes 

In order to study the limiting null distribution of the test statistics S g , n in (12.61) and T g ^ n 
in (12.71) . it is necessary to investigate the asymptotics of sequential empirical processes 
constructed from PWM estimators. The main result of this section is used in the proof, 
given in Appendix [Dl of Proposition 12.31 but it might also be of independent interest. 

As we consider a slightly more general theoretical framework than in the main sections 
of the paper, a richer notation is necessary. Let X be a random variable with c.d.f. F 
and let iq,..., u p be p functions on [0,1]. The reals f3 v = E\ Xu{F{X) lh v G {u i, ..., v p }, 
assuming they exist, are called generalized PWM following Diebolt et al. (2008). When 
the functions zy are defined as Vi(x) = x ri (l — x) Si , x G [0,1], for some integers r* > 0 
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and Si > 0, one recovers a very usefull subset of “classical” PWM that can be used 


to fit a large number of d i stribu tions (see, e.g., lOreenwood et all Il979t iHosking et al 


19851: IHosking and WallisL I1987T) . In the sequel, we write u = (zq,..., v p ) and /3„ = 


(Pvi ) • • • J ) • 


Assuming that we have at hand a random sample X\,..., X n from an unknown 
c.d.f. F, a natural nonparametric estimator of j3 u based on a subsample X^,..., Xi, 
1 < k < l < n, is 

/3u,k:l • j fii/p,k\l)i (C.l) 


where f3 Vu k-.i stands for j3 i}k -.i dehned in ( 12 .3p with the difference that 7 in ( 12 . 2 p is set to 
zero. 


For any x > 0, let |ycj be the largest integer smaller or equal than x. Furthermore, 
recall that A = {(s,£) G [0, l] 2 : s < t) and that A n (s,t) = {\nt\ — \ns\)/n, ( s,t ) G A. 
The aim of this section is to establish the weak convergence of the vector of empirical 

processes M„ jn = (M^,..., M^J, where, for any v G {v 1 ,..., v p }, 


(s, t) \frl\ n {^S, f) [fiv, |nsj+l:|ntj (®> ^ A. 


(C.2) 


The desired result is established under the following condition. 

Condition C.l. The functions zq,. .., v p are continuously differentiable on [0,1]. Fur¬ 
thermore, Xi ,..., X n are i.i.d. from a continuous distribution with c.d.f. F such that: 

(i) fi v . = E[X 1 u i {F(Xi)}] < 00 for all i G {1,..., p}, 

(n) E([A 1 ^{F(X 1 )} - /3„J [X\Vj{F{Xf)} - &.]) < 00 for all i,j G {1,... ,p}, 

(Hi) there exists a G [0,1/2) such that sup^gg \H a (x)\ < 00 , where 

H a (x) = x[F(x){l — F(x)}] a , iGl. 


It is easy to verify that if one of the functions zq,..., v p is constant on [0,1], then (i) 
and (ii) in Condition 1C. II are equivalent to requiring that E(Xf) < 00 . 

Proposition C.2 (Limiting null distribution of Under Condition \C. 11 


hdj , n I • • • , J-VJljyp jj 


'W M,. = 


LL^l , . . . , -LVJi^ p 


in {£°°(A)} P , where M„ lin ,..., are defined in (1C.21) and M^,..., M Fp 
Gaussian processes with continuous trajectories such that 


are centered 


vi(s,t) =M„.(0 ,t) — My.(0, s), (s,t) G A, j G {1 ,...,p}, 


(C.3) 


and 


Cov{M^.(0, s), M^.(0, t)} = (s At)Cov(Y u .,Y Uj ), ( s,t ) G A, i,j G {l,...,p}, (C.4) 

where, for any 1 / G {zq,..., u p }, Y v is defined in (12.9H . 
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Proof. Fix v G {ui,, u p } and consider the decomposition 


M Uj n(s, t) = X^ n (s, t ) + Y U)U (s, t ), 


(M) G A, 


where 


\/nX n (s, t) (/3i/,[nsj+l: [ntj fiv, [_nsj+l: pitj ) > 

Y u,n(.Sp) ^~TlX n (s, t) (/3i/ 5 [nsj+l:[>iij Pv)i 


(C.5) 


with 


fiv, [nsj-t-npitj 


[ntj — [nsj 


\nt\ 

Y. X tV {F(X t )}. 


(C.6) 


i=[nsj+l 

In addition, let U\,... ,U n be the unobservable sample obtained from Ah,, X n by the 
probability integral transformations Ui = F^Xf), i G {1,... , n}, and recall the definition 
of the empirical process Wq,^ in (IB.211 and its weak limit W a given in Proposition IB.11 
The first step of the proof is to show that 


(Y„ 1>n ,..., Y„ n , W a, n ) (Y^,..., Y v , W a ) 


(C.7) 


in {^°°(A)}p x £°°(A x [0,1]). 

The proof of the convergence in distribution of the finite-dimensional distributions of 
(Y„ ljn ,..., Y Up)1l , W a t n) to those of (Y^,..., Y Up , W a ) is omitted for the sake of brevity. 
To prove (1C.7ft . it remains thus to show joint asymptotic tightness, which is implied 
by marginal asymptotic tightness. From Proposition IB.II we have that W a<n W Q in 
7°°(A x [0,1]). Next, fix v G { v\ ,..., u p } and notice that, since of = Var[Ad^'{i ? (Ad)}] < 
oo, we have, from the functional central limit theorem, that Y„ jn (0, •) Y^(0, •) in 

^°°([0,1]), where Y„(0, •) = oyU with U a standard Brownian motion. Since, for any 
(s, t) G A, Y Vt n(s, t) = Y^ n (0, t) — Y I/)Tl (0, s), an immediate consequence of the continuous 
mapping theorem is that Y^ n -w Y„ in ^°°(A), where, for any (s,t) G A, Y v (s,t) = 
oy{U(t) — U(s)}. This completes the proof of (1C. 7ft . 

Fix again v G {iq,..., v p }. Since sup xeR \H a (x)\ < oo, sup ue r 01 i \v'(u)\ < oo and F is 
continuous, the map : 7°°(A x [0,1]) —> £°°(A) defined as 


■, a {f)(s,t)= / F{F(x)}H a (x)f{s,t,F(x)}dF(x), (s,t) e A, / G T(A x [0,1]), 

Jr 

(C.8) 

is continuous. The weak convergence in (1C.7K and the continuous mapping theorem then 
imply that 


Yj 1 ,n+^,a(Wa, I! ), • • • , '^v p ,n,-\-'ipv p ,a (W ain )^ ^Y Vl +^ liQ (W Q ), . . . , Y t , p +^ p>a (W Q ) 


in {£°°(A)} P . To complete the proof, it remains thus to show that, for any v G {zq,..., u p }, 


sup \X Ujn (s,t) -^„, a (W a)n )(s,i)}| = op(l), (C.9) 

(s,t)e A 
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where X„ )n is defined in (1C.51) . Fix v G { 17 ,..., z/ p } and consider the empirical process 
X' (s,t) = / i/{F(a;)}//' a (x)W airi {s,t,F(x)}dF Lns j + i :Lnt j(x), (s,t)eA, (C. 10 ) 


where TjnsJ+ 1: lMl is defined in (12.2[) with 7 = 0 . To prove (1C.9[) . we shall first show that 

(C.ll) 


and then that 


snp |X„ )tl (s,t) -X' (s,t)| = o p (1), 

(s,t)e A 


SU P | X *,n(M) -^,a(W ain )(s,t)| = 0 P (1). 
(s,t)e A 


(C.12) 


Proof of (1C.111) . On one hand, starting from (1C.51) . we have that, for any (s,t) G A, 

^ IMI 

X UtU (s,t) = —7= ^ X i [^{F\ns\ + l-.\nt\{Xi)} - u{F(Xi)}] 

i=[ns\-\-l 

= / xy/n\ n (s,t) [^{F LnsJ+ i. LntJ (a;)} - u{F(x)}] dF LtMj+ i : L„ t j(a:). 

Jr 

On the other hand, (IC.lOp can be rewritten as 

X L,*»(M) = / ^ , {F 1 (x)}v / nA n (s,t){F LnsJ+ i :Lnt j(a:) - F(x)}dF LnsJ+ i :Ln ij(a;), (s,t)eA. 

Jr 

Then, the snpremnm on the left of (1C. Ill) is smaller than 

sup y/nX n (s,t) / |x| \v{F [nsi+1:lnti (x)} - v{F(x)} 

(s,t)eA Jr 

—v'{F(x)}{F[ns\+l:lnt\ (x) ~ -F(x)}| dF|_ res j+i : \_ nt \ {x). (C.13) 
By the mean value theorem, for any iGR, there exists U* stx G (0,1) such that 

v{ F \na\+l:lnt\ &)} ~ v{F(x)} = t/[F{x) + U* stx {F\ |n S J+l:|ntJ (x) ~ ^0)}] 

x {-F|nsJ+ 1: M ( x ) 

Hence, (1C.131) is smaller than 


sup / |z| I u'[F(x) + U* tS ^ x {F lns j+i :Lnt j(x) - F(x)}} - u'{F(x)}\ 

(s,t)£ A J M 


X | \Zri\ n {^S) [nt\ (*^) |cl-^^LnsJ+1:L n ^J (*^) ’ 


that is, 


sup / \H a (x) \ | v'[F(x) + U* jSjtjX {F lnsi+1:[nti (x) - ^(X)}] - ^{F(x)}\ 
(s,t) gA Jr 


X |W ai „{s, t, F(x) }|dF LnaJ+1: L„ t j (®) • 
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Therefore, 


sup \X^ n (s,t) -X^ n (s,t)| 

(s,t)e A 

< sup |t/'[F(x) + C/* Siti;r {F Lns j +1:Lnt j(a;)-F(x)}]-^{F(a:)}| 

(s,£,:e)gAxR 

x sup |W Q)n {s, £, F(a;)}| x sup \H a (x)\. (C.14) 

(s,t,:r)G Axl fcGR 

Fix e, r/ > 0. From the previous inequality, since sup ug [ 0jl ] |z/(w)| < oo, sup xgR \H a (x)\ < 
oo, and using the fact that W a , n (s,t,-) = 0 when s = t and the asymptotic uniform 
equicontinuity in probability of W a>n (a consequence of Proposition IB.IF there exists 
5 G (0,1) such that, for all sufficiently large n, 


P sup IX I/ ,„(s,t) - X' un (s,t)\ > £ 

(s,t)e A 
t—s<8 


< P ( 2 sup |z/(it)| x sup \H a (x)\ x sup |W ain (s, t, u)\ > e ) < rj/2. 

\ ue[0,l] xeR (s,t,u)eAx[o,i] 

To prove (1C.lip . it remains therefore to show that, for all sufficiently large n, 


P sup |X vn (s,t) -Xj'(s,t)j > £ < r]/2. 


(s,t)e a 
t—s>5 


Let A s = {(s, t) G A : t — s > <5}. To show the latter, it suffices to show that 


The fact that 


sup |X I/ , n (s,t) -X' l/n (s,f)| = o P (l). 

(s,()eA J 

in £°°(A x [0,1]), where B n is defined in (IB.IF implies that 


sup |G|n S j + i:|ni|(u) “ u \ = SU P |i^LnsJ+l:(ntj OO “ F(x) \ = Op (1.), (C.15) 

(s,t,w)GA 5 x [0,1] (s,t,:r)G A 5 xR 

since G Lris j + i : L nt j{F(x)} = (x) for all (s,t,x) gAxB. The previous display 

then implies that 

sup | v'[F(x) + U*^ x {F lnsi+Hnti {x) - F(x)}] - v'{F(x)}\ = o P (l) 

(s,t,x) GA 5 xR 


since 1 / is (uniformly) continuous on [0,1]. The latter combined with (1C.141) in which the 
suprema are restricted to (s,t) G A 5 and the fact that sup( Sit)X \ eAxR |W Q)n {s, t, F(x)}| = 
0 P (1) implies finally (1C. lip . 

Proof of (1C.121) . We have 

SUp \X, n ( S ’t) ~ ^MW a , n )(M)}| < SUp 1X^(5, f)| + Slip |^, a (W Q , n )(s,t)| 

(s,t)G A (s,t)eA (s,t)£A 

<2 sup |i/(?z)| x sup \H a (x)\ x sup |W ajn (s, t, n)|. 

ue[0,l] XSR (s,t,u)gAx[0,l] 
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Fix e,r] > 0. Proceeding as for the proof of (1C.121) . from the previous inequality, the 
fact that W a ,n{s,t,-) = 0 when s = t and the asymptotic uniform equicontinuity in 
probability of W a ,m there exists 5 G (0,1) such that, for all sufficiently large n, 


P sup |X[, n (s,£)-^,a(W Q)n )(s,t)}| > £ <7]/2 , 

t—s<5 


and thus, to prove (1C. 121) . it remains to show that 


sup |X' (s,t) - ^,a(W Q)n )(s,t)| = Op(l), 

(s,t)eA { 


(C.16) 


where A 6 = {(s,t) G A : t — s > 5}. Proving the previous display is equivalent to showing 
that X^ n — ^,a(W a , n ) converges weakly to (s,t) i->- 0 in ^(A 5 ). Fix (s,t) G A 3 and let 
G(u) = u for all u G [0,1]. From (IC.15j) . (W Q)n (s, t, •), G^j+nintj) -w (W a (s, t, •), G) in 
{£°°([0, l])} 2 , which implies that 


x ^ v'{F(x)}H a (x)W atn {s,t, F(x)}, F [nsi+1:[nti 

x ^ u'{F(x)}H a (x)W a {s, t, F(x)}, F 


in {^(R)} 2 , where R = [—00,00]. The assumptions of the proposition ensure that 
x 1 —y u'{F(x)}H a (x) is a continuous function on R and therefore that the process 


x 1— y i/ / {F(a;)}i/ Q ,(a;)W Q ,{s, t, F(x)} G £° 


has co ntinuous trajectories. Combining the latter fact with Lemma 3 of (Holmes et ah 
( 2013 ) and the continuity of the map V-Va defined in (1C.81) . we obtain, from the continuous 
mapping theorem, that X' un (s,t) — Vv,a(W ai „)(s, t) 0. The latter weak convergence 
being equivalent to convergence in probability to zero, we have that the finite-dimensional 
distributions of X(, n — ^(W^n) converge weakly to those of (s,t) 1 —> 0 in £°°(A 5 ). 

It remains to show asymptotic tightness. To do so, we shall prove the || • ||i-asymptotic 
uniform equicontinuity in probability of X'„ n — ^ a (W Qi „) in £°°(A 5 ). By the triangle 
inequality, it suffices to show the latter for X], n only as converges weakly 

to ijj„ n fW n ) in £°°(A S ), the limit having continuous trajectories. By Problem 2.1.5 in 
van der Vaart and Wellnerl (120001 ). it remains thus to prove that, for any positive sequence 

-i- 0, 

sup |x;„(s, i) - x;„(s', 01 = op(l). (C.17) 

(s,t),(s / ,£ / )G 

\s — s'| + |£— t'\<Q n 

The supremum on the left of the previous display is smaller than A n + B n , where 


A n = sup 


u'{F(x)}H a (x)W a>n {s, t, F(x)} d-Fjnsj +i: |_ ntJ (x) 

- [ F{F(x)}H a (x)W ain {s',t' : F(x)}dF lnsi+1:lnti (x) 
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and 


B n = sup 

(s,t),(s / ,t / )G 
\s — s'| + |t — t'\<an 


v'{F(x)}H a (x)W a , n {s', t', F(x)}dF lnsi+1:[nt] (x) 

- [ F{F(x)}H a (x)W ajn {s , ,t , ,F(x)}dF lns/i+1:[nt/i (x) 


Concerning A n , we have 

A n < sup \i/(u) \ x sup \H q (x)\ x sup |W ai „(s, t, u) — W Q)n (s / , t', u)\ = op(l). 

uG[0,l] zER (s,t),(s' ,t')eA s ,ue[o,i] 

|s—s'l + lt— t'\<a n 

The term B„ is smaller than 


sup 

\s — s'l + lt—t'|<a n 


[nt] 


\nt\ — [nsj 


Lnf'J 


u’{F(X i )}H a {X^^ n {s' 1 t' 1 F(X i )} 


i=[nsj+l 


+ 


sup 

(s,i),(s / ,£ / )G 
\s — s'l + li — i'|<a n 


- ^ I /{F(X i )}^ a (X i )W a , n {5 , ,£ / ,F(X i )} 

i=[7is'j+l 

\nt'\ 

Y v'{F( y X i )}H a (X i )'W ol}n {s', t', F(Xj)} 


[nt\ — [ns] [nt’\ — l ns 'l 


i=\ns' J+l 

which is smaller than 

I \. n t\ — \nt'\ | + | |nsj — [ns'J | 


2 x sup 

(s,t),(s',t')e a 5 
\s — s'l + lt — t'lKan 


\nt\ — [nsj 


x sup |z/(u)| 
«e[o,i] 


x sup \H a (x)\ x sup |W Q)n (s, £, u)\ = Op(1). 

zeR (s,t, u)eAx[0,l] 


Hence, (1C. 171) and thus (1C. 161) hold, which implies (IC.12|) . 

Finally, (IC.3j) is a consequence of the fact that, for any v G {zq,..., u p }, 


(s,t) =Y„(s,t) + / u'{F(x)}H a (x)W a {s,t, F(x)}dF(x) 


and the analogue properties for Y v and W a , while (1C.41) follows from standard calcula¬ 
tions. ■ 


D Proof of Proposition 12.3 

Proposition 12.31 is a corollary of a slightly more general result that we shall first prove 
using Proposition 1C. 21 shown in Appendix O To be consistent with Appendix [Cj we keep 
using the richer notation introduced therein. Let zq,..., v p be p functions on [0,1], let 
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g be a function from M p to M defined on D containing (3„ = (/3 Ul ,... ,/3 Up ), and consider 
the generic statistic 


Sv g n = max 

l<k<n -1 


k(n — k) 1 
n 3 / 2 


(A/,l:fc £ P, is,k-\-l:n ^ ^(/^i^rfe) di.ftiSyk+l.n') 


where /3^1-fc and /3„,fc+i :ri are defined analogously to (1C. ID . When p = 3, ui(x) = 1, 
z/ 2 (a;) = x, z/ 3 (x) = x 2 , x G [0,1], P = in (ES]) and g G {s^, fitr, fi^}, is 

in (12H1) . 

Proposition D.l (Limiting null distribution of *S' I y i9 , n ). If Condition \C. 1\ holds. E(|Ad|) < 
oo, P is an open convex subset ofW p containing (3 U and g is continuously differentiable 
onP, S Vi g in converges in distribution to a Uj9 sup s£ [ 01 ] |U(s)—sU(l)|, where U is a standard 
Brownian motion on [0,1] and 

v 

°l, g = Y1 d i9((3u)djg{^)Coy(Y Ui ,Y Vj ), 
i,j =1 


with Y u . defined in (12.9p . 

Proof. For any s G [0,1], let 

®i',g,n('S) \/^A n (0, s) A n (s, 1) 1 (/^^p |^j s J G P, @v,\ns\+l:n ^ P) ftv,l:\ns\) 9 iftu, [nsj +l:n) }> 

and notice that 

Su, g ,n= sup |D Vlffl n(s)|. (D.l) 

se[o,i] 

Furthermore, let 

fdliv,g,n(S) t) \fn\ n (yS, t) {g(ftv,\ns\+l-.\nt \) ^(/^i/)}, (s, t) G A. 

Under H 0 in (11.21) implied by Condition 1C. 11 it is easy to verify that 

D>iy,g,)i(s) = A n (s, l)l(/3 Iy> i : |_ ns j G P, ftu,\ns\+l'.n £ I l )M|/ 1 g,ti(0, s) 

-^(0,8)10 G P)M„g }Tl (s, 1), sG [0,1]. (D.2) 

To show the desired result, we shall first prove that 


A. n SUp 1 [ftv,l-.\ns\ G P, A/,|nsJ+l:n G P) 

seto,i] 

and that 


B n SUp Ifyftu,l:\ns\ G P, ftu,\ns\+l\n eP) 
se[o,i] 


/,S,n(0, s) - dig(Pv) M I/il „(0, s) 


2=1 


s,g : n(Sj 1) - ^^(/3 I/ )M t , i>n (s, 1) 


2=1 


where 


^1,72 5 * * * 5 lvlL Z/p 


in are defined in (IC.2p . 


- o P (l) 
(D.3) 

= Op(l), 
(D.4) 
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Let us start with (11 ).•'>[) . Using the mean value theorem when the indicators are equal 
to one, A n is smaller than A ijU , where 

A,n = sup l(/3„,i :Lns j G V)\dig(f3* 1:[nsJ ) - 0i0(/3„)||M„. in (O, s)|, (D.5) 

se[o,i] 


and where /3* 1 .. , lies on the segment between (3„ and i:\ns\- To show (II). -'ll) , it suffices 

therefore to prove that A itn = Op(l) for all i G {1,... ,p}. 

Fix v G {zq,..., u p }. Since j3 u < oo, from the strong law of large number, we have 
that ——>■ fiv, where [3 V) i, n is defined analogously to (1C.6|) . Thus, 


Pu,l:n ~ Pu,l:n < SUp \u{F l:n (x)} ~ u{F(x)}\ X - V] |A^| 0 

a;gR n ^ 


%— 1 


since E(|Ad|) < oo, and the supremum on the right converges to 0 almost surely as a 
consequence of the Glivenko-Cantelli lemma and the continuous mapping theorem. It 
follows that $ u ,i:n /3 U for all u G {zq,..., u p }, and therefore that (3 v ^, n — 

Fix i G {1,... ,p}. The continuous mapping theorem next implies that dig(/3^^ :n ) 
dig((3 u ), which in turn implies that sup fc>1 G T , )|<9jg(/3 iyi i : / c )— dig((3 u )\ < oo almost 

surely. 

Let g > 0. Then, there exists M > 0 such that 


fc>i 


P Z sup 1 i/Vi:).. e ©) - d,g(/3 r ) > M \ < rj/4. 


(D.6) 


Furthermore, let e > 0 and consider e/M > 0. Using the asymptotic uniform equiconti- 
nuity in probability of (a consequence of Proposition 1C.2ft . there exists 5 G (0,1) 

such that, for all sufficiently large n, 


P < sup 

se[o,<5) 


<,n(0.s)l > e/M }■ < r]/4. 


(D.7) 


Then, 


P (Ai, n > e) < P < sup 1(^,1;^] G V)\dig(/3* 1:[nsi ) - dig{p v )\\M v .^, s)| > e 
1 se[o,5) 

+ P<{ sup 1(^,1: L»wj G V)\dig(PZ , 1:L n S j) - ^(/3^)l|Mv i ,n(0,s)|} > £ . (D.8) 

The first term on the right is smaller, for all sufficiently large n, than 


P <{ sup 1 e V) dig((3 Vtl;k ) - dig((3 v ) 

k> 1 


> M \ +P { sup 

se[o,5) 


i(0, s)| > e/M ^ < rj/2 , 
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where the last inequality follows by (ID. 6(1 and (ID. 7(1 . The second term on the right 
of (ID. 81) is smaller, for all sufficiently large n, than rj/2. Indeed, 


sup l(p„ M ns} e V)\dig(Pl 1:M ) - diS(/3„)||M I/ . in (0,s)| 
*£[<5,1] 


< sup \dig0t Mnsi ) - dig(f3„)\ x sup |M Vij „(0, s)| 
*£[<5>1] *6[0,1] 


Op(l) 


since the second supremum is Op(l) from the weak convergence of M y . )n , while the hrst 
supremum is op(l) as a consequence of the fact that, for any i G {1,... ,p}, 

sup \Ki, 1:LH - P"i\ =ri~ 1/2 sup A n (0, s) _1 x sup |M Vi)n (0, s)| = o P (l) 

se[<5,l] *£[<5,1] sG[0,1] 

and the continuous mapping theorem. Hence, A ijn defined in (ID . 5|) converges in proba¬ 
bility to zero for all i G {1,... ,p}, which implies that (ID .3[) holds. The fact that (ID. 41) 
holds follows immediately since B n in (ID. 41) (written as a maximum) is nothing else 
than A n in (ID. 31) computed from the sequence X n ,..., Xi. Therefore, for any e > 0, 
P (B n > e) = P (A n > e) —> 0 by (ID. 31) . 

Next, let 

^u,g,n = S Up G T), /3i/,[ns\+l:n *7 l®i;,g,ra( S ) 15 

*£[0,1] 

where 

p 

D ^ 9 ,n( s ) = 5^^(/3 I/ ){A n (s,l)M I/ . jn (0,s) - A n (0,s)M I/ . in (s, 1)}, s G [0,1]. (D.9) 

i =1 


Then, starting from (ID.II) and (ID.21) . and using the reverse triangle inequality, we obtain 
that 


Ik _ Q' 

\ u i',g,n u u,g,n\ 


< A n + B n — Op(l), 
where A n and B n are defined in (ID.3j) and (ID.41) . 


(D.10) 


We shall now verify that 
where 


Q' 

nr/r 


-^ |B = op(l), 

= sup |B' (s)|. 

*£[ 0 , 1 ] 


( 0 . 11 ) 

( 0 . 12 ) 


To do so, we use the decomposition S' u gn — S'” gn 


q/ _ qn _i_ qn _ qm 


where 


q// 


sup l(/3„,[ns|+l:n ^ X>)|D^ (s)|. 

*£[ 0 , 1 ] 


Since (3 u ^-.k —A /3„ as k —> oo, we have that, for almost every u, there exists k u such 
that k>k u implies that /3„ ; i :fc (u;) G V. Hence, for almost every u, 


S 'v,g,n(u) ~ S” g n (uj) = sup l{/3„ lL „ s j + i:nM G V}\B)' v g n ( S , u)\ 

*£[fc^/n, 1] 


Slip l{A,,Lmj+l:nM € V}\W v ^ n (s,u)\ -> 0, 
*£[ 0 , 1 ] 
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which implies that S' ug n — S'" = op(l). Similarly, we obtain that 

SUp G V)\W (s) \ ~ S'” = 0 P (1). 

se[o,i] 


Then, using the fact that S'^ gn - S'” gn is sup se[0>1] l(A,,i: Lns j e v )\ D L, P ,n 
computed from the sequence X n ,..., Xi, we obtain that, for every e > 0, 


(a) I 


n/f/ 


n\s:, g , n - s" g j >e) = p 


sup 1(A,,1: L „ sJ g V)\®'„ gn (s)\ - S'" 
se[0,i] 


> £ > ->■ 0. 


Therefore, (ID.Ill) holds. Next, from Proposition 1C.21 and the continuous mapping theo¬ 
rem, we have that D ' vgn W v g in f°°([0,1]), where D ' vgn is defined in flD.Ojl . 

v P 


i= 1 


2—1 


for s G [0,1], and M V1 ,..., Mj, are defined in Proposition lC.2l The desired convergence in 
distribution of S v>g!n finally follows from (ID.lOp . (ID.111) . (ID . 12[) and standard calculations. 


Proof of Proposition 12.31 The convergence in distribution of S 9tn is an immediate 
consequence of Proposition ID.II It remains to prove that S g ^ n — = op(l). The latter 

is a consequence of the fact that, for any i G {1, 2, 3}, 

In = sup IM^OjS) - V / ^A n (0,s)(6i ! i : Lnq - AJ | = o P (l) (D.13) 

se[o,i] 

and 

Jn SUp 1) ■\ZnA n (s, 1) (^i,[nsj+l:n AO I Op(l)> (D.14) 

se[o,i] 

where 0,1: [nsj and 0,|nsj+i :n are defined analogously to (12.41) . The above statements are 
clearly true for i — 1 as bi^i = Ai,fc : z for all integers 1 < k, l < n. For z = 2, it can be 
verihed that the supremum on the left of (ID . 13(1 restricted to s G [2/n, 1] can be rewritten 
as 


sup —= 

sG[2/n,l] v ^ 


[nsj 


L?isJ 


J2 X i F l:lns}( X i) ~J2 X i 


|7TsJ F\. 1 


i =1 


i= 1 


[nsj — 1 

[nsj 


< 


sup 


^ sG[2/n,l] [ris] 


—7 T |A r i| < —p sup jPlXi 

- 1 fe ^ k fef 


since E(|Ab|) < oo. The case i — 3 is similar. To show (ID.14(1 . we use the fact that J n is 
I n computed from the sequence X n ,... ,X 1 and conclude from (ID . 131) . 

Starting from (ID. 13)1 and (ID. 14ft , it is possible to adapt the proof of Proposition ID. II 
to show that S g>n — T g ^ n = o P (l). We omit the details for the sake of brevity. ■ 
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Table 1: Percentage of rejection of Hq computed from 1000 samples of size n e 
{50,100,200,400} generated from the GEV(0,1,£). 


PWM / T 9)n GPWM / S h , n PWM / T g , n GPWM / S h , n 


e 

n 

k 

9 a 

k 

hf,, 

h(j 

h 

e 

n 

k 

9 a 

k 

hg , 

ha 

h 

- 1.0 

50 

3.0 

0.7 

4.1 

7.2 

0.8 

0.2 

0.1 

50 

3.8 

3.2 

4.0 

1.2 

0.6 

1.4 


100 

2.9 

1.6 

3.7 

7.4 

1.7 

0.2 


100 

4.9 

3.4 

2.7 

2.7 

1.4 

1.3 


200 

4.0 

2.9 

4.1 

5.4 

3.2 

0.6 


200 

4.7 

4.2 

2.5 

3.8 

2.4 

2.9 


400 

3.8 

3.3 

4.4 

3.9 

2.9 

1.1 


400 

4.4 

4.5 

3.4 

4.1 

3.5 

3.9 

- 0.9 

50 

2.6 

0.7 

3.8 

4.8 

0.2 

0.4 

0.2 

50 

6.2 

5.4 

4.4 

2.3 

1.5 

1.5 


100 

4.3 

2.0 

4.2 

6.3 

1.5 

0.2 


100 

5.2 

4.9 

3.6 

3.6 

2.2 

1.6 


200 

4.8 

2.8 

5.1 

4.8 

2.4 

0.7 


200 

5.0 

4.4 

3.2 

4.3 

2.6 

4.2 


400 

4.0 

3.3 

3.6 

3.4 

3.5 

1.1 


400 

4.6 

4.2 

4.2 

4.2 

3.9 

4.4 

- 0.8 

50 

2.2 

0.8 

4.0 

4.5 

0.2 

0.3 

0.3 

50 

5.7 

3.7 

5.0 

2.2 

0.8 

1.4 


100 

2.2 

2.4 

2.8 

4.3 

1.6 

0.1 


100 

4.3 

4.1 

2.9 

2.6 

1.4 

2.0 


200 

4.1 

3.8 

3.8 

4.7 

3.6 

0.5 


200 

6.1 

6.0 

2.5 

4.5 

3.0 

3.2 


400 

5.6 

4.1 

4.2 

5.9 

4.3 

2.0 


400 

5.0 

4.8 

2.6 

3.7 

3.7 

4.0 

- 0.7 

50 

2.4 

1.4 

3.9 

3.6 

0.1 

0.0 

0.4 

50 

6.6 

4.4 

8.1 

3.1 

1.5 

2.1 


100 

3.1 

2.5 

4.3 

4.2 

1.1 

0.1 


100 

6.8 

5.6 

4.5 

4.5 

2.0 

2.9 


200 

3.9 

3.7 

4.1 

5.0 

2.7 

0.6 


200 

6.1 

5.0 

3.3 

4.0 

3.8 

3.0 


400 

3.7 

3.7 

3.9 

4.3 

3.5 

1.3 


400 

5.3 

4.5 

2.1 

3.6 

4.8 

4.3 

- 0.6 

50 

2.2 

1.9 

3.6 

2.2 

0.1 

0.3 

0.5 

50 

8.0 

4.2 

9.5 

2.7 

1.3 

3.5 


100 

4.1 

2.9 

3.8 

3.7 

0.7 

0.3 


100 

9.2 

7.0 

6.8 

4.8 

3.2 

2.9 


200 

5.1 

4.1 

3.5 

4.3 

2.5 

0.7 


200 

5.8 

4.5 

3.4 

4.8 

3.3 

3.5 


400 

4.3 

4.6 

4.0 

4.4 

4.1 

2.4 


400 

5.4 

5.6 

2.9 

4.5 

4.1 

4.2 

- 0.5 

50 

3.3 

3.2 

4.5 

2.7 

0.3 

0.7 

0.6 

50 

6.7 

4.8 

16.6 

4.2 

2.5 

3.4 


100 

3.4 

2.3 

4.7 

3.7 

0.5 

0.7 


100 

8.5 

6.8 

10.9 

5.6 

4.4 

3.0 


200 

3.4 

3.8 

4.9 

3.8 

2.2 

1.3 


200 

6.4 

5.9 

7.6 

4.5 

3.9 

3.7 


400 

4.5 

4.4 

4.1 

4.5 

3.3 

2.4 


400 

5.3 

7.6 

2.7 

3.8 

4.8 

3.4 

- 0.4 

50 

2.5 

3.2 

4.9 

2.0 

0.4 

0.3 

0.7 

50 

7.8 

4.0 

18.5 

3.9 

2.1 

3.2 


100 

3.6 

3.5 

3.9 

2.8 

1.2 

1.0 


100 

10.4 

7.2 

17.1 

5.8 

5.1 

2.4 


200 

3.0 

5.0 

3.9 

3.2 

2.7 

1.5 


200 

7.1 

6.6 

8.5 

5.3 

4.7 

3.2 


400 

5.3 

4.9 

3.7 

5.2 

3.5 

1.9 


400 

7.4 

6.9 

6.5 

6.3 

5.4 

3.5 

- 0.3 

50 

2.8 

4.2 

4.4 

1.3 

0.7 

0.6 

0.8 

50 

9.6 

3.1 

27.4 

4.6 

4.3 

5.9 


100 

4.6 

3.0 

4.3 

3.9 

0.9 

0.7 


100 

10.1 

7.3 

24.3 

7.2 

6.0 

3.2 


200 

4.6 

4.2 

4.6 

4.4 

1.9 

2.1 


200 

9.4 

7.7 

18.8 

6.8 

6.3 

2.7 


400 

4.2 

4.1 

4.7 

3.6 

3.0 

3.1 


400 

7.5 

7.4 

10.6 

5.4 

5.1 

3.0 

- 0.2 

50 

2.6 

2.9 

3.0 

1.5 

0.4 

0.4 

0.9 

50 

11.9 

5.5 

31.1 

7.0 

5.3 

6.1 


100 

5.1 

4.3 

3.9 

3.6 

0.9 

1.0 


100 

10.8 

7.4 

33.6 

7.8 

6.0 

2.2 


200 

4.7 

4.6 

4.3 

4.0 

2.7 

2.0 


200 

10.4 

9.3 

25.5 

7.4 

6.0 

3.0 


400 

5.1 

4.1 

4.6 

3.9 

3.0 

3.2 


400 

7.6 

7.4 

18.9 

6.0 

6.7 

2.7 

- 0.1 

50 

4.0 

3.8 

3.9 

1.8 

0.7 

0.8 

1.0 

50 

12.0 

4.6 

37.8 

7.6 

6.4 

9.5 


100 

4.2 

4.2 

4.7 

2.6 

1.6 

1.5 


100 

11.6 

7.4 

39.1 

7.4 

6.7 

4.2 


200 

4.7 

4.4 

4.5 

4.1 

2.4 

2.1 


200 

12.3 

8.5 

32.4 

9.6 

10.1 

3.0 


400 

4.0 

3.7 

4.2 

3.4 

2.4 

3.8 


400 

7.8 

5.2 

27.5 

5.8 

6.2 

2.8 

0.0 

50 

4.7 

3.0 

4.0 

1.8 

0.3 

1.3 

1.1 

50 

13.9 

4.8 

38.8 

10.8 

7.6 

12.2 


100 

3.6 

3.3 

2.8 

2.6 

1.3 

1.5 


100 

15.3 

8.3 

44.9 

12.0 

9.6 

5.3 


200 

4.5 

3.6 

3.8 

3.6 

2.2 

1.7 


200 

9.8 

7.4 

44.4 

7.9 

9.5 

2.3 


400 

4.6 

3.4 

3.7 

4.0 

2.8 

4.0 


400 

7.5 

6.4 

38.0 

7.0 

8.3 

2.4 
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Table 2: Percentage of rejection of Hq computed from 1000 samples of n E 
{50,100,200,400} independent block maxima obtained from blocks of size 1 from the 
GPD(0,1,Q._ 


£ 

n 

PWM / T 9)U 

GPWM / S h>n 

£ 

n 

PWM / T g , n 

GPWM / S h , n 

k 

9 a 

k 


Hq- 

h 

k 

9(7 

k 

V 

h(j 

h 

- 1.0 

50 

3.0 

3.9 

1.7 

1.6 

3.0 

1.8 

0.1 

50 

6.1 

2.9 

5.5 

3.5 

1.5 

5.0 


100 

3.7 

4.3 

2.5 

2.2 

2.8 

1.6 


100 

5.7 

3.3 

4.3 

3.9 

2.2 

4.2 


200 

3.8 

4.3 

3.1 

2.7 

3.6 

2.9 


200 

4.9 

3.3 

2.4 

3.8 

2.9 

3.5 


400 

4.2 

4.7 

3.7 

3.9 

4.3 

3.1 


400 

5.5 

3.5 

4.2 

4.8 

3.7 

4.9 

- 0.9 

50 

3.4 

3.5 

2.8 

1.5 

3.0 

1.8 

0.2 

50 

6.6 

3.4 

5.0 

2.6 

1.5 

3.9 


100 

4.9 

3.5 

3.8 

2.8 

1.8 

2.5 


100 

5.4 

3.6 

3.7 

3.6 

2.4 

4.3 


200 

4.4 

4.4 

3.5 

3.6 

3.7 

3.2 


200 

5.7 

4.0 

3.1 

5.0 

4.5 

4.8 


400 

3.5 

4.3 

3.5 

3.3 

3.6 

4.2 


400 

5.0 

5.3 

3.1 

3.8 

4.6 

4.4 

- 0.8 

50 

3.8 

1.6 

2.3 

1.1 

1.3 

1.7 

0.3 

50 

6.3 

3.7 

6.3 

3.0 

1.8 

5.2 


100 

3.5 

2.7 

3.0 

2.7 

1.9 

2.0 


100 

6.5 

4.8 

5.5 

4.4 

2.9 

4.6 


200 

4.4 

3.7 

3.5 

3.5 

2.9 

3.4 


200 

5.7 

5.0 

2.3 

4.0 

3.4 

4.7 


400 

4.8 

4.2 

4.3 

4.4 

3.5 

4.5 


400 

5.2 

4.1 

2.8 

4.3 

3.8 

4.5 

- 0.7 

50 

4.5 

1.9 

2.0 

1.9 

1.0 

1.8 

0.4 

50 

9.8 

5.1 

9.0 

5.7 

2.9 

4.3 


100 

5.2 

2.5 

3.2 

3.5 

2.0 

1.5 


100 

7.8 

6.1 

5.0 

5.2 

3.0 

4.1 


200 

4.3 

3.3 

3.8 

2.9 

2.1 

3.8 


200 

7.3 

6.0 

3.3 

6.1 

5.0 

3.6 


400 

3.4 

4.4 

5.3 

3.4 

3.9 

4.6 


400 

5.9 

4.3 

2.8 

5.4 

3.9 

4.4 

- 0.6 

50 

4.1 

2.2 

2.5 

2.0 

1.8 

2.4 

0.5 

50 

7.8 

3.9 

11.0 

4.3 

2.0 

3.8 


100 

4.5 

3.1 

3.3 

3.5 

2.2 

3.8 


100 

7.0 

6.8 

9.2 

5.1 

4.9 

4.5 


200 

4.0 

3.2 

3.3 

2.8 

2.9 

2.9 


200 

6.8 

5.0 

6.2 

4.8 

3.1 

4.3 


400 

4.1 

4.5 

3.4 

4.1 

3.8 

4.0 


400 

6.5 

5.9 

3.5 

4.3 

4.6 

3.3 

- 0.5 

50 

4.6 

1.9 

2.3 

2.1 

1.2 

1.8 

0.6 

50 

8.4 

4.5 

17.0 

5.1 

2.8 

5.5 


100 

4.1 

2.7 

2.4 

2.3 

1.9 

2.5 


100 

10.1 

7.1 

11.1 

6.7 

4.9 

4.9 


200 

6.3 

3.6 

5.4 

4.5 

3.6 

3.7 


200 

9.3 

7.6 

7.0 

6.2 

5.3 

3.6 


400 

5.7 

2.9 

4.9 

5.0 

3.0 

5.6 


400 

5.9 

6.4 

4.8 

5.1 

5.3 

3.8 

- 0.4 

50 

5.0 

2.1 

3.7 

1.9 

1.0 

3.1 

0.7 

50 

9.9 

5.7 

22.1 

6.2 

3.9 

6.7 


100 

4.5 

1.7 

3.5 

2.7 

1.5 

3.1 


100 

9.9 

7.0 

19.8 

7.2 

5.2 

3.5 


200 

4.6 

3.3 

3.3 

3.7 

2.9 

3.4 


200 

8.0 

7.5 

11.2 

6.6 

5.4 

3.2 


400 

5.2 

3.8 

3.8 

4.6 

3.8 

3.3 


400 

7.1 

5.5 

5.2 

6.5 

4.4 

4.4 

- 0.3 

50 

3.3 

1.9 

3.6 

1.3 

0.6 

1.8 

0.8 

50 

11.4 

5.9 

27.4 

7.1 

5.5 

6.9 


100 

5.1 

3.4 

3.2 

2.9 

2.2 

3.2 


100 

13.2 

7.1 

25.2 

8.8 

6.6 

3.8 


200 

5.4 

2.7 

3.0 

3.8 

2.8 

4.0 


200 

10.4 

7.2 

16.7 

7.7 

7.5 

2.5 


400 

6.0 

3.6 

4.1 

5.8 

3.4 

4.1 


400 

6.3 

6.2 

12.0 

5.1 

4.8 

3.2 

- 0.2 

50 

4.9 

1.4 

3.5 

1.8 

0.8 

2.1 

0.9 

50 

12.4 

5.9 

31.0 

9.9 

7.7 

9.4 


100 

4.5 

2.6 

3.3 

3.0 

1.3 

2.8 


100 

12.4 

8.6 

32.0 

7.9 

8.9 

5.2 


200 

5.2 

3.1 

5.0 

4.0 

2.7 

4.1 


200 

10.2 

6.8 

25.0 

8.1 

7.2 

1.7 


400 

4.2 

3.9 

3.2 

4.7 

4.1 

4.0 


400 

8.1 

7.0 

18.5 

6.3 

5.6 

2.4 

- 0.1 

50 

6.8 

2.8 

3.8 

4.0 

0.9 

3.5 

1.0 

50 

12.6 

4.9 

37.0 

8.5 

7.0 

9.1 


100 

5.1 

2.2 

3.5 

3.3 

1.8 

4.0 


100 

13.6 

6.8 

40.4 

11.1 

8.8 

4.8 


200 

4.6 

3.7 

3.0 

4.2 

3.0 

3.8 


200 

10.9 

7.4 

35.1 

9.7 

9.5 

2.1 


400 

5.5 

4.6 

5.5 

5.3 

3.9 

4.8 


400 

8.0 

7.4 

25.6 

6.7 

7.3 

3.1 

0.0 

50 

6.1 

3.5 

2.9 

2.9 

1.8 

3.7 

1.1 

50 

16.3 

5.6 

40.0 

11.7 

9.2 

14.3 


100 

5.0 

3.1 

4.4 

3.6 

1.5 

4.2 


100 

15.6 

8.3 

47.3 

10.8 

11.3 

5.7 


200 

6.2 

3.7 

4.3 

5.1 

4.0 

4.1 


200 

11.2 

8.3 

41.5 

11.0 

9.6 

3.7 


400 

5.7 

4.7 

4.5 

5.1 

4.0 

5.1 


400 

8.7 

7.7 

37.1 

8.0 

7.4 

2.5 
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Table 3: Percentage of rejection of H 0 computed from 1000 samples of size n e {100, 200} 
such that the \nt\ first observations are from a GEV(0,l,-0.4) and the n — \nt\ last 
observations are from a GEV(0,1,£). 

PWM / T 9)U GPWM / S h>n 


£ 

n 

t 

F n 


s 2 

d n 

k 

9<j 

k 

k 

fl(j 

k 

-0.2 

100 

0.25 

5.7 

5.1 

4.3 

3.2 

3.8 

10.2 

2.3 

1.7 

1.2 



0.75 

6.1 

4.9 

7.0 

4.5 

3.4 

5.1 

3.3 

0.9 

1.4 



0.50 

5.5 

6.5 

8.1 

4.2 

3.4 

12.5 

2.9 

1.7 

2.0 


200 

0.25 

5.2 

7.0 

6.0 

3.7 

4.0 

17.7 

3.7 

2.2 

5.1 



0.75 

7.3 

8.3 

8.5 

5.5 

3.9 

15.3 

4.5 

1.9 

3.5 



0.50 

8.5 

11.7 

10.7 

4.9 

3.5 

29.4 

4.7 

1.6 

6.9 

0.0 

100 

0.25 

5.8 

7.6 

4.5 

4.8 

4.2 

23.1 

2.6 

1.2 

3.2 



0.75 

7.6 

12.6 

13.0 

5.3 

4.4 

11.6 

4.2 

1.1 

4.5 



0.50 

7.8 

16.8 

10.2 

4.4 

2.6 

32.8 

3.0 

0.8 

9.1 


200 

0.25 

8.3 

16.8 

7.6 

4.5 

2.6 

54.4 

3.8 

1.1 

22.7 



0.75 

8.0 

22.2 

29.8 

5.7 

5.2 

40.2 

3.7 

2.1 

15.4 



0.50 

10.3 

32.1 

30.7 

4.6 

3.1 

76.1 

3.2 

1.5 

37.9 

0.2 

100 

0.25 

8.0 

9.7 

2.4 

6.1 

4.7 

42.7 

3.1 

1.4 

14.3 



0.75 

8.3 

23.6 

11.7 

5.8 

5.5 

16.0 

4.2 

1.5 

9.2 



0.50 

12.0 

27.7 

10.3 

6.8 

6.7 

52.8 

3.8 

3.0 

26.8 


200 

0.25 

12.0 

24.9 

4.9 

7.1 

4.7 

78.2 

5.0 

2.8 

58.2 



0.75 

11.5 

51.3 

39.4 

6.8 

4.4 

59.7 

4.2 

2.2 

44.6 



0.50 

22.0 

65.9 

35.0 

6.6 

6.4 

92.4 

4.9 

3.2 

81.8 

0.4 

100 

0.25 

9.6 

11.1 

0.5 

9.6 

7.3 

62.6 

4.3 

2.0 

27.5 



0.75 

9.3 

34.4 

11.7 

6.4 

9.5 

21.7 

3.5 

2.5 

20.9 



0.50 

15.7 

40.8 

7.2 

8.1 

10.1 

69.5 

4.2 

2.5 

50.3 


200 

0.25 

15.5 

25.1 

2.3 

11.1 

7.0 

89.2 

4.6 

2.9 

85.1 



0.75 

13.7 

71.4 

36.3 

7.9 

7.6 

63.9 

3.3 

4.2 

79.2 



0.50 

29.6 

80.6 

26.9 

12.3 

8.7 

95.7 

6.3 

3.4 

97.4 
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Table 4: Percentage of rejection of H$ computed from 1000 samples of size n e {100, 200} 
such that the \_nt\ first observations are from a GEV(0,1,0) and the n — \ nt\ last obser¬ 
vations are from a GEV(0,1,£). 

PWM / T 9)n GPWM / S h>n 


£ 

n 

t 

F n 


s 2 

9 /-t 

9<j 

k 

V 

h(j 

k 

0.2 

100 

0.25 

5.7 

4.4 

1.2 

5.8 

4.2 

6.7 

2.8 

2.0 

4.1 



0.75 

6.7 

7.4 

3.2 

4.6 

4.6 

3.9 

3.5 

2.2 

2.9 



0.50 

4.1 

5.5 

2.2 

4.6 

4.8 

6.7 

2.9 

1.3 

4.3 


200 

0.25 

8.2 

6.7 

2.2 

5.9 

4.4 

10.0 

4.8 

2.4 

8.8 



0.75 

5.7 

9.9 

7.8 

5.7 

4.8 

6.0 

4.4 

3.3 

7.2 



0.50 

6.1 

9.9 

7.3 

5.2 

4.8 

14.7 

4.6 

2.8 

15.9 

0.4 

100 

0.25 

5.9 

4.6 

0.4 

6.2 

4.4 

15.5 

3.3 

2.1 

7.0 



0.75 

6.1 

11.3 

2.5 

7.0 

8.0 

5.9 

3.8 

1.7 

6.6 



0.50 

6.4 

13.1 

4.4 

6.7 

8.0 

19.2 

2.9 

2.2 

17.4 


200 

0.25 

7.1 

8.0 

1.8 

5.4 

5.2 

28.8 

3.9 

2.4 

27.9 



0.75 

7.7 

25.9 

15.2 

6.1 

7.2 

13.0 

4.5 

3.2 

24.5 



0.50 

7.6 

30.6 

10.3 

5.1 

6.4 

37.9 

3.7 

2.8 

54.2 

0.6 

100 

0.25 

6.7 

5.1 

0.7 

8.6 

7.9 

39.9 

4.5 

3.2 

17.8 



0.75 

7.3 

19.5 

5.2 

7.5 

13.7 

13.1 

3.9 

3.6 

15.0 



0.50 

7.7 

17.4 

1.6 

7.4 

12.5 

36.8 

3.8 

3.4 

35.1 


200 

0.25 

8.4 

9.4 

1.3 

9.1 

8.8 

57.6 

5.2 

3.9 

57.7 



0.75 

8.3 

47.4 

21.9 

7.1 

15.0 

24.3 

4.0 

5.6 

56.6 



0.50 

15.4 

45.7 

11.5 

10.6 

14.7 

68.3 

5.1 

5.7 

84.1 
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Table 5 : Percentage of rejection of H 0 computed from 1000 samples of size n e { 100 , 200 } 
such that the \ nt \ first observations are from a GEV ( 0 , 1 , 0 . 2 ) and the n — \ nt \ last 
observations are from a GEV ( 0 , 1 ,£). 

PWM / T ff ,„ GPWM / S h>n 


e 

n 

t 

F n 


s 2 

9 n 

9 cr 

k 

hfj , 

Hey 

k 

0.4 

100 

0.25 

5.7 

2.8 

0.5 

6.8 

4.7 

6.1 

4.0 

1.2 

4.5 



0.75 

6.1 

4.3 

1.1 

6.1 

5.5 

3.7 

4.3 

2.3 

3.6 



0.50 

5.9 

6.5 

0.3 

6.7 

6.0 

6.6 

3.8 

1.7 

5.0 


200 

0.25 

5.2 

4.4 

0.8 

5.5 

4.8 

7.8 

4.0 

3.2 

8.9 



0.75 

6.0 

7.9 

4.5 

5.9 

6.5 

3.4 

4.8 

4.1 

8.7 



0.50 

5.4 

8.9 

2.9 

4.3 

5.0 

9.7 

3.2 

2.3 

15.9 

0.6 

100 

0.25 

6.7 

2.4 

0.2 

8.1 

6.3 

22.6 

4.5 

2.9 

10.7 



0.75 

7.3 

8.1 

2.3 

8.1 

9.4 

7.0 

5.3 

3.0 

6.7 



0.50 

7.6 

7.8 

1.3 

9.7 

9.4 

18.4 

4.6 

3.5 

16.3 


200 

0.25 

6.4 

6.1 

0.5 

6.9 

7.6 

27.7 

4.8 

4.6 

24.5 



0.75 

7.0 

21.0 

7.4 

6.7 

10.4 

9.4 

4.5 

5.7 

23.5 



0.50 

7.5 

19.2 

4.5 

7.7 

10.0 

29.5 

4.0 

3.7 

48.9 

0.8 

100 

0.25 

6.0 

3.1 

0.3 

10.7 

8.4 

49.5 

4.9 

4.1 

16.5 



0.75 

5.8 

13.9 

2.5 

8.8 

17.0 

13.6 

4.4 

4.1 

11.1 



0.50 

8.3 

13.1 

1.0 

12.9 

18.4 

43.9 

5.3 

5.8 

28.3 


200 

0.25 

9.2 

5.6 

0.8 

13.4 

11.3 

58.6 

6.5 

5.6 

46.4 



0.75 

7.4 

38.3 

11.2 

9.7 

23.0 

23.6 

4.6 

7.1 

49.1 



0.50 

14.2 

30.6 

4.0 

15.8 

24.7 

60.5 

6.9 

8.5 

75.3 
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Table 6: Percentage of rejection of Hq computed from 1000 samples of size n e {100, 200} 
such that the [nt\ first observations are from a GEV(0,0.5,£) and the n — \nt\ last 
observations are from a GEV(0,1,£). 

PWM / T g>n GPWM / S h>n 


£ 

n 

t 

F n 


Sn 

9n 

9cr 

k 

hn 

ha 

h 

- 0.4 

100 

0.25 

17.0 

7.3 

54.3 

4.5 

62.2 

9.8 

4.8 

38.1 

1.8 



0.75 

20.1 

13.4 

94.5 

13.5 

82.3 

2.6 

13.9 

58.9 

0.9 



0.50 

40.7 

12.4 

98.1 

8.5 

93.6 

5.7 

9.8 

81.9 

1.0 


200 

0.25 

50.0 

8.4 

97.8 

5.5 

97.4 

11.8 

4.2 

93.7 

3.5 



0.75 

44.7 

19.6 

99.9 

15.5 

98.7 

2.8 

15.6 

96.5 

0.4 



0.50 

85.6 

24.4 

100.0 

10.5 

100.0 

5.1 

10.9 

100.0 

0.9 

0.0 

100 

0.25 

15.0 

9.0 

18.3 

6.8 

57.6 

6.7 

5.8 

39.1 

2.1 



0.75 

13.7 

17.7 

59.6 

12.5 

78.6 

3.1 

10.3 

50.6 

0.7 



0.50 

26.0 

20.7 

64.7 

11.0 

91.7 

5.5 

9.9 

79.1 

0.5 


200 

0.25 

28.1 

15.7 

50.6 

6.3 

95.6 

7.5 

4.8 

89.9 

2.7 



0.75 

29.5 

33.1 

91.0 

15.0 

97.9 

1.7 

13.9 

95.8 

0.7 



0.50 

63.1 

41.6 

94.6 

10.8 

99.9 

4.9 

10.5 

99.8 

0.7 

0.4 

100 

0.25 

8.2 

3.9 

0.6 

6.9 

25.2 

7.1 

5.0 

22.0 

2.6 



0.75 

10.9 

11.0 

5.2 

14.3 

58.1 

6.4 

12.6 

42.3 

1.6 



0.50 

17.6 

14.5 

5.9 

12.2 

70.3 

7.9 

10.7 

63.6 

1.6 


200 

0.25 

18.8 

9.9 

2.7 

7.3 

62.4 

6.4 

6.0 

69.9 

5.4 



0.75 

20.4 

22.6 

11.8 

12.9 

85.9 

3.6 

13.2 

85.5 

1.3 



0.50 

47.5 

29.7 

12.2 

11.3 

94.5 

6.1 

10.7 

95.2 

2.2 

0.8 

100 

0.25 

7.3 

2.7 

0.0 

9.4 

11.3 

25.6 

6.6 

14.1 

3.9 



0.75 

10.1 

5.0 

0.6 

18.9 

34.1 

23.4 

14.0 

36.1 

2.9 



0.50 

12.7 

4.4 

0.3 

13.1 

32.6 

23.0 

10.4 

45.8 

2.8 


200 

0.25 

16.8 

3.0 

0.6 

8.7 

23.8 

21.6 

6.7 

40.4 

3.6 



0.75 

13.2 

7.8 

1.0 

17.3 

55.8 

14.6 

14.9 

70.6 

1.2 



0.50 

36.8 

9.5 

0.9 

15.7 

61.5 

18.6 

13.4 

79.9 

2.1 
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Table 7: Percentage of rejection of Hq computed from 1000 samples of size n G {100, 200} 
such that the \_ nt \ first observations are from a GEV(0,1,£) and the n — [ nt \ last obser¬ 
vations are from a GEV(0.5,1,£). 

PWM / T g>n GPWM / 5 ft , n 



n 

t 

F n 

•E n 

Sn 

9 n 

9 cj 

k 

h/i 

fl(j 

h 

- 0.4 

100 

0.25 

39.3 

39.2 

2.2 

27.4 

3.9 

3.9 

27.6 

1.1 

0.5 



0.75 

36.5 

38.4 

3.9 

27.9 

3.4 

3.8 

19.3 

1.0 

0.5 



0.50 

58.0 

61.9 

2.1 

48.0 

3.5 

4.1 

42.7 

0.6 

0.5 


200 

0.25 

67.6 

72.1 

3.6 

57.7 

3.2 

4.5 

53.9 

1.6 

1.5 



0.75 

68.6 

73.3 

4.1 

57.1 

5.1 

4.2 

47.8 

2.5 

1.0 



0.50 

88.6 

91.0 

1.2 

79.3 

3.0 

3.7 

74.1 

1.4 

2.0 

0.0 

100 

0.25 

28.2 

19.7 

3.7 

26.5 

3.4 

2.6 

23.8 

0.9 

1.1 



0.75 

28.2 

23.3 

2.2 

29.0 

4.7 

4.2 

21.2 

1.8 

2.0 



0.50 

48.6 

39.0 

1.9 

48.3 

2.3 

3.4 

42.1 

0.7 

1.7 


200 

0.25 

52.8 

42.2 

4.0 

54.6 

4.1 

2.5 

52.8 

2.3 

1.8 



0.75 

51.7 

42.4 

1.7 

54.9 

4.5 

4.0 

48.1 

3.6 

2.7 



0.50 

76.7 

66.6 

3.5 

78.1 

3.7 

3.1 

76.1 

2.4 

2.3 

0.4 

100 

0.25 

23.1 

4.7 

0.5 

22.1 

4.7 

4.7 

19.8 

2.0 

2.3 



0.75 

22.8 

7.8 

0.5 

27.2 

5.2 

4.3 

22.0 

2.8 

3.1 



0.50 

43.5 

10.6 

0.3 

45.2 

3.8 

3.6 

42.0 

1.4 

3.0 


200 

0.25 

48.8 

8.5 

1.2 

48.4 

5.0 

3.2 

51.9 

3.0 

2.7 



0.75 

48.6 

8.5 

1.0 

54.1 

5.6 

2.8 

49.9 

3.9 

4.6 



0.50 

74.0 

16.8 

0.9 

71.6 

4.7 

2.3 

76.6 

2.8 

3.1 

0.8 

100 

0.25 

22.9 

1.2 

0.2 

21.7 

7.1 

21.0 

18.5 

4.9 

2.5 



0.75 

24.2 

2.5 

0.2 

27.9 

6.7 

19.7 

26.5 

5.3 

4.0 



0.50 

46.2 

2.9 

0.2 

38.9 

6.8 

23.2 

41.1 

4.9 

1.3 


200 

0.25 

48.6 

2.0 

0.4 

40.0 

7.9 

15.4 

45.6 

6.9 

2.2 



0.75 

52.0 

3.3 

0.4 

46.4 

7.9 

17.5 

51.8 

6.9 

3.9 



0.50 

77.5 

3.4 

0.4 

64.7 

7.8 

13.6 

71.6 

5.1 

2.0 
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Tabic 8: Minimum and maximum of <7 M (&i :n ), 9a(bi-. n ) and g^(bi :n ) as well as minimum and maximum of the p-values of the tests 
based on T g ^ n in (12.7ft with g respectively equal to g^, g a , g^ for 1000 different de-tied samples generated from each dataset as 
explained in Section [5j The first two columns provide the names of the R packages containing the datasets and the names of the 
latter. The columns n and n' report the number of block maxima and the number of distinct block maxima, respectively. 


Package 

Dataset 

n 

n' 



Estimates 



p-values of the test based on T g ^ n with 


9cr(Pv,l:n) 

9i (/3i/,l:n) 

9 = 


9 = 

: 9cr 

9 = 

-k 

min 

max 

min 

max 

min 

max 

min 

max 

min 

max 

min 

max 

evd 

lisbon 

30 

21 

95.79 

96.22 

12.62 

13.07 

-0.16 

-0.13 

0.152 

0.205 

0.167 

0.271 

0.416 

0.630 

evd 

oxford 

80 

19 

84.23 

84.46 

4.20 

4.44 

-0.34 

-0.26 

0.099 

0.248 

0.534 

1.000 

0.413 

1.000 

extRemes 

HEAT (Tmax) 

43 

9 

112.89 

113.22 

1.91 

2.32 

-0.41 

-0.21 

0.002 

0.029 

0.590 

1.000 

0.766 

1.000 

extRemes 

HEAT (-Tmin) 

43 

15 

-70.70 

-70.37 

3.67 

4.07 

-0.18 

-0.09 

0.000 

0.002 

0.130 

0.529 

0.173 

0.898 

extRemes 

ftcanmax 

100 

78 

135.75 

135.96 

55.51 

55.79 

0.13 

0.13 

0.724 

0.757 

0.430 

0.467 

1.000 

1.000 

ismev 

fremantle 

86 

31 

1.49 

1.49 

0.14 

0.14 

-0.21 

-0.19 

0.006 

0.009 

0.433 

0.578 

1.000 

1.000 

ismev 

portpirie 

65 

42 

3.88 

3.88 

0.20 

0.20 

-0.06 

-0.04 

0.537 

0.603 

0.788 

0.949 

0.782 

0.928 

ClusterMax 

Station 8 (St-Lizier) 

228 

58 

1.97 

1.98 

1.95 

1.97 

0.14 

0.15 

0.217 

0.261 

0.296 

0.378 

0.068 

0.091 

ClusterMax 

Station 14 (Aurillac) 

228 

73 

2.39 

2.40 

2.87 

2.89 

0.29 

0.29 

0.152 

0.174 

0.622 

0.683 

0.124 

0.138 

ClusterMax 

Station 57 (Nevers) 

228 

56 

1.82 

1.84 

1.90 

1.92 

0.18 

0.18 

0.287 

0.349 

0.155 

0.202 

0.034 

0.046 

ClusterMax 

Station 78 (Niort) 

228 

59 

1.96 

1.97 

2.15 

2.17 

0.22 

0.22 

0.038 

0.047 

0.016 

0.022 

0.448 

0.523 













